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PREFACE 


This  report  provides  documentation  of  three  computer  programs  for 
performing  settlement  analysis  of  foundations  and  embankments.  The 
report  was  written  as  part  of  the  normal  operation  of  the  joint  U.  S. 

Army  Engineer  Waterways  Experiment  Station  (WES)  and  U.  S.  Army  Engineer 
Division,  Lower  Mississippi  Valley,  Computer  Center  for  Fiscal  Years 
1978  and  1979. 

The  three  computer  programs  documented  herein  are  10016,  MAGSETII 
and  FD31.  Program  10016,  which  was  developed  by  Mr.  Douglas  Spaulding, 
Foundation,  Materials,  and  Survey  Branch,  St.  Paul  District,  determines 
vertical  stresses  beneath  footings  and  embankments.  MAGSETII  was 
written  by  Messrs.  R.  L.  Schiffman,  D.  M.  Jubenville,  and  V.  Partvka  of 
the  University  of  Colorado  to  calculate  the  magnitudes  of  settlement  of 
multilayered  soil  systems.  Dr.  Roy  E.  Olson,  University  of  Texas,  Austin, 
developed  FD31  to  determine  time-settlement  relationships  for  cohesive 
soils  due  to  large  uniformly  distributed  loads. 

The  documentation  was  put  together  in  a  package,  with  example  runs 
and  comparisons  with  hand  computations,  by  Mr.  Reed  L.  Mosher,  Computer- 
Aided  Design  Group,  Automatic  Data  Processing  (ADP)  Center,  WES,  under 
the  direct  supervision  of  Dr.  N.  Radhakr ishnan.  Special  Technical  Assis¬ 
tant,  ADP  Center.  This  report  was  written  by  Mr.  Mosher  and  Dr.  Radha- 
krishnan.  Mr.  D.  L.  Neumann  was  Chief  of  the  ADP  Center. 

COL  J.  L.  Cannon,  CE,  and  COL  N.  P.  Conover,  CE,  were  Directors  of 
WES  during  the  preparation  of  this  report.  Mr.  F.  R.  Brown  was  Technical 
Director. 
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CONVERSION  FACTORS,  INCH-POUND  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


Inch-pound  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 


To  Obtain 


feet  0.3048 

kips  (1000  lb  force)  4.448222 

kips  (force)  per  square  foot  47.88026 
pounds  (force)  per  square  foot  47.88026 
pounds  (mass)  per  cubic  foot  16.01846 
square  feet  0.09290304 


metres 

kilonewtons 

kilopascals 

pascals 

kilograms  per  cubic  metre 
square  metres 
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PROGRAM  INFORMATION 


This  settlement  package  described  herein  is  operational  on  the  U.  S. 

Army  Engineer  Waterways  Experiment  Station's  Honeywell  0-635  time¬ 
sharing  system  at  Vicksburg,  Miss.,  and  on  the  Office  of  Personnel  Man¬ 
agement's  Honeywell  6000  Series  time-sharing  system  at  Macon,  Oa.  The 
file  names  used  for  the  programs  are  listed  below  with  short  descriptions 
of  how  to  access  each.  It  is  assumed  that  the  user  knows  how  to  sign  on 
to  the  system  he  is  using. 

10016 

*  FORT 

*  RUN  WESLIB/CORPS/ 10016, R 

MAGSETII 

*  FORT 

*  RUN  WESLIB/CORPS/ 10010, R 

FD31 

*  FORT 

*  RUN  WESLIB/CORPS/ 10011 ,R 


COMPUTER  PROGRAMS  FOR  SETTLEMENT  ANALYSIS 


PART  I:  INTRODUCTION 


Purpose 


1.  A  package  of  three  programs  for  performing  settlement  analysis 
of  foundations  and  embankments  is  documented  in  this  report.  The  pro¬ 
grams  are  based  on  theories  and  methods  accepted  by  practicing  engineers 
and  presented  in  universities  throughout  the  United  States. 

2.  The  package  can  be  a  very  powerful  and  time  saving  aid  to  the 
foundation  engineer  in  the  analysis  of  complex  foundation  systems.  With¬ 
out  the  use  of  the  computer,  solutions  to  problems  involving  such  systems 
could  be  lengthy  and  tedious  and  leave  room  for  human  error.  The  pro¬ 
grams  allow  the  foundation  engineer  to  be  more  creative  by  providing 
more  time  to  explore  innovative  alternatives. 


Programs  in  the  Package 


3.  The  package  consists  of  three  separata  programs:  10016, 
MAGSETII,  and  FD31.  10016  determines  vertical  stresses  beneath  footings 

and  embankments.  It  was  developed  by  Douglas  Spaulding  (1968)  of  the 
St.  Paul  District.  MAGSETII  calculates  the  magnitudes  of  settlement 
of  multilayered  soil  systems.  It  was  written  by  R.  L.  Schiffman,  D.  M. 
Jubenville,  and  V.  Partyka  (1976)  at  the  University  of  Colorado.  Addi¬ 
tions  to  the  program  to  compute  the  degree  of  consolidation  have  been 
made.  FD31  develops  time-settlement  relationships  for  cohesive  soils 
due  to  large  uniformly  distributed  loads.  It  was  written  by  Roy  Olson 
at  the  University  of  Texas  at  Austin. 


A.  This  report  provides  documentation  for  the  three  computer 
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programs  used  in  the  package.  Theories  on  which  the  programs  are  based, 
along  with  capabilities  of  the  programs,  are  discussed.  Input /output 
for  the  programs  is  discussed  using  three  example  problems.  One  of  the 
example  problems  is  taken  from  Engineer  Manual  1110-2-1904  (Headquarters, 
Department  of  the  Army,  Office  of  the  Chief  of  Engineers  1953).  Docu¬ 
mentation  for  the  programs,  as  provided  by  the  original  authors,  is  ref¬ 
erenced  in  this  report.  Original  documentation  for  the  programs  can  be 
obtained  from  the  Engineering  Computer  Program  Library  (ECPL)  at  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES). 


PART  II:  METHODS  AND  CAPABILITIES 


5.  Program  10016  can  calculate  vertical  stress  distributions  in 
a  soil  profile  based  on  either  Boussinesq  or  Westergaard  solutions. 

Both  methods  assume  that  the  soil  is  homogeneous,  isotropic,  and  linearly 
elastic.  Westergaard  further  assumes  that  there  are  no  lateral  deforma¬ 
tions.  These  assumptions  do  not  completely  model  actual  soil  behavior, 
but  without  these  assumptions  solutions  are  only  possible  using  more 
sophisticated  techniques.  In  most  cases,  the  results  obtained  using 
these  simplified  assumptions  are  reasonably  accurate  when  compared  to 
field  observations  (see  Spaulding  1968).  10016  allows  the  user  to 

analyze  rectangular  loadings  and/or  embankment  loadings  in  a  three- 
dimensional  layout.  The  embankment  loadings  are  applied  by  number  of 
uniform  rectangular  shaped  layers  with  the  width  decreasing  with  height. 
This  allows  the  user  to  consider  problems  involving  time-dependent  loads 
due  to  construction,  etc.  The  user  has  the  option  to  choose  the  horizon¬ 
tal  or  vertical  plane  to  be  investigated.  The  capabilities  are  illus¬ 
trated  best  in  the  example  problems  presented  in  Part  III. 


Program  MAGSETII 


6.  Program  MAGSETII  utilizes  Terzaghi's  one-dimensional  consoli¬ 
dation  theory,  simplified  to  apply  to  a  two-dimensional  condition,  for 
estimating  settlements  in  cohesive  soils.  The  effective  stress  history 
for  each  layer  or  for  the  total  profile  can  be  input  to  the  program. 

The  program  applies  a  vertical  stress  influence  factor,  due  to  the  load¬ 
ing,  to  the  effective  stress  history.  Under  this  effective  stress  his¬ 
tory,  some  very  complex  loadings  can  be  accounted  for,  such  as:  unload¬ 
ing  due  to  excavation,  temporary  and/or  permanent  changes  in  water  table, 
live  loads  applied  to  the  structure,  and  loadings  due  to  adjacent  struc¬ 
tures  or  construction.  Granular  soils  are  handled  by  empirical  correla¬ 
tions  to  static  or  dynamic  penetration  field  tests.  MAGSETII  takes  into 
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account  strain  influence  with  depth  in  granular  soils.  It  has  two 
built-in  methods  to  account  for  strain  influence  (Figure  1),  or  the 
user  can  enter  a  set  of  influence  factors.  Also,  for  granular  soils, 
three  methods  are  available  for  estimating  settlements:  Meyerhof's, 
D'Appolonia' s,  and  Schmertmann ' s.  The  first  two  methods  use  data  from 
a  standard  penetration  test;  Schmertmann ' s  method  uses  data  from  a  static 
cone  penetrometer  test  (see  Schuffman,  Jubenville,  and  Partyka  1976). 


Figure  1.  Strain  influence  in  sands 


7.  A  subroutine  to  compute  rate  of  settlement  has  been  added. 
These  computations  are  based  on  Terzaghi's  theory  and  methods  described 
in  EM  1110-2-1904.  This  addition  gives  the  user  the  option  to  consider 
the  effect  of  time-dependent  loading  on  the  rate  of  settlement,  as 
outlined  in  EM  1110-2-1904. 


Program  FD31 

8.  Settlement  and  time-settlement  relationships  for  compressible 
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materials  are  computed  in  FD31  based  on  Terzaghi's  one-dimensional  con¬ 
solidation  theory.  The  program  is  only  valid  for  one-dimensional  analy¬ 
sis.  The  differential  equations  derived  from  the  theory  are  solved  by 
finite  difference  methods  of  analysis.  (See  Olson.)  FD31  provides  the 
user  with  a  very  versatile  tool  to  compute  settlement  and  rate  of  settle¬ 
ment  for  cohesive  soil.  The  program  allows  for  a  stratified  soil  pro¬ 
file,  subject  to  time-dependent  loadings;  soils  that  are  not  linearly 
elastic,  which  may  be  subject  to  large  and  nonuniform  strains,  and  whose 
coefficients  of  permeability  and  compressibility  may  vary  with  effective 
stresses;  and  stress  conditions  that  are  altered  by  a  changing  water 
table  and  settlement-dependent  submergence  of  the  soil. 

9.  FD31  is  a  specialized  program.  It  is  very  sensitive  to  the 
data  input,  and  the  user  must  be  careful  in  correctly  modeling  the  in 
situ  situation.  Input  data  come  from  standard  laboratory  tests  and 
field  observations.  The  program  does  not  take  into  account  the  influence 
of  vertical  stress  distribution  with  depth. 

Loads 


10.  Geometric  configurations  play  an  important  part  in  the  choice 
of  program.  Two  basic  types  of  loadings  can  be  handled.  These  are: 

(a)  concentrated  loads  and  (b)  uniformly  distributed  loads.  If  the 
width  of  the  structure  applying  the  load  at  the  surface  is  relatively 
small  in  comparison  to  the  depth  of  the  compressible  soil,  the  load  can 
be  considered  to  be  concentrated.  Loading  conditions  which  fall  under 
this  category  are:  strip  footings,  spread  footings,  some  raft  founda¬ 
tions,  and  also  embankments  in  which  the  base  is  relatively  small  in 
comparison  to  the  compressible  soil  being  considered.  If  the  area  being 
loaded  is  wide  compared  to  the  depth  of  the  compressible  soil,  the  load 
should  be  considered  as  uniformly  distributed.  Loading  conditions  which 
fall  under  this  category  are:  fills,  embankments,  and  large  excavations. 

11.  For  analysis  of  concentrated  loads,  MAGSETII  is  the  primary 
program  used.  It  can  handle  a  multilayered  soil  profile  of  cohesive 
and/or  granular  material.  It  can  account  for  preloadings  and  unloadings. 
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When  estimating  settlements  for  cohesive  material  under  a  concentrated 
loading,  10016  is  used  to  calculate  the  vertical  stress  distribution 
beneath  the  point  being  investigated.  The  data  from  this  program  can 
ba  used  directly  in  MAGSETII.  To  achieve  the  best  accuracy,  large  layers 
of  compressible  material  are  subdivided  into  several  smaller  layers. 

12.  For  analysis  of  large  uniformly  distributed  loads  on  layers 
of  compressible  material,  FD31  is  used.  In  the  case  where  a  compressible 
soil  and  a  granular  soil  are  in  the  same  profile,  the  settlement  due  to 
the  granular  material  would  be  negligible  in  comparison  with  the  settle¬ 
ment  of  the  cohesive  soil. 
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PART  III:  EXAMPLE  PROBLEMS  ILLUSTRATING  INPUT/OUTPUT  FOR 
PROGRAMS  10016  AND  MAGSETII 


13.  In  this  Part,  two  example  problems  are  solved  using  programs 
10016  and  MAGSETII.  Input/output  for  the  two  programs  is  also  described. 
Results  of  problem  1  are  compared  with  hand  solutions.  Problem  2  is 
taken  from  EM  1110-2-1904,  and  the  results  are  compared  with  values  from 
that  source. 


Example  Problem  1 

Program  10010 

14.  Organization  of  input.  The  input  data  are  categorized  into 
three  groups:  header  lines,  loading  configuration  data,  and  stress  dis¬ 
tribution  data.  The  first  of  these  groups  consists  of  five  lines  of  data 
describing  the  particular  run.  The  second  group  describes  the  geometric 
configuration  and  loads  applied  by  embankments  and/or  footings.  The 
third  group  defines  the  type  of  analysis  (Boussinesq  or  Westergaard) 

and  the  location  and  direction  (whether  distribution  along  a  vertical  or 
horizontal  plane  is  desired) .  The  amount  of  data  required  for  the 
second  and  third  data  groups  is  dependent  on  the  complexity  of  the 
problem  and  the  output  required. 

15.  Mode  of  input.  Input  to  the  program  can  be  either  from  the 
terminal  or  from  a  data  file.  (Example  problem  1  was  solved  using  data 
input  from  the  terminal;  problem  2,  which  is  discussed  later  in  this 
Part,  was  solved  by  reading  data  already  stored  in  a  data  file.)  All 
input  is  in  free  field.  Data  items  can  be  separated  by  a  blank  space 
or  a  comma.  If  information  is  being  read  from  a  data  file,  each  line 
of  data  must  be  preceded  by  a  line  number.  When  operating  from  the 
terminal,  the  program  can  create  files  to  save  the  input  data  and  the 
output.  Detailed  input  with  definitions  of  input  variables  for  program 
10016  is  shown  in  later  paragraphs  of  this  Part  using  problem  1  as  an 
example . 

16.  Problem  definition.  Figure  2  shows  a  plan  view  of  two  rectan¬ 
gular  footings  loading  the  soil  profile  shown  in  Figure  3.  The  profile 
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Figure  2.  Plan  view  of  footings 

consists  of  10  ft*  of  fill  material  and  15-,  6.5-,  and  20.5-ft  layers 
of  clay  material  with  a  25-ft  layer  of  sand  and  gravel  sandwiched  be¬ 
tween  the  last  two  clay  layers.  The  water  table  is  25  ft  below  the 
surface. 

17.  The  footings  are  placed  after  excavating  10  ft  of  material. 

Then  10  ft  of  new  fill  material  is  placed  and  compacted.  As  construc- 

2 

tion  continues,  the  structure  applies  loads  of  2.0  and  2.5  kips/ft  , 
respectively,  to  the  footings.  At  the  end  of  construction,  0.5  kip/ft'' 
is  relieved  from  the  footings.  Table  1  shows  this  information  and  the 
times  these  events  occur. 

18.  Figures  4-6  show  void  ratio  versus  effective  stress  curves 


*  A  table  of  factors  for  converting  inch-pound  units  of  measurement 
to  metric  (SI)  units  is  presented  on  page  3. 
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Table  1 

Loading  Conditions  for  Example  Problem  1 


Load  Increment 
1 
2 

3 

4 

5 


_ Load _ 

10  ft  of  fill  excavated 

10  ft  of  new  fill 
2 

2-kips/ft  loading 
2 

2.5-kips/ft  loading 
2 

0.5-kip/ft  unloading 


Time  Interval 
days _ 

0  to  50 

50  to  75 

75  to  200 

200  to  300 

300  to  350 


for  the  three  clay  layers  in  the  soil  profile.  This  information  is 
given  in  tabular  form  in  Table  2.  Table  3  shows  results  from  standard 
penetration  tests  for  the  sand  layer. 

19.  For  this  example,  the  settlement  is  estimated  under  the  center 
of  the  footing.  Program  10016  is  first  used  to  calculate  the  vertical 
stress  influence  factors  beneath  the  center  of  the  footing  for  the  cohe¬ 
sive  layers.  This  will  be  done  at  the  midheight  of  each  layer. 

20.  Input.  Data  required  for  program  10016,  arranged  by  groups, 
are  shown  below. 

a.  Problem  information.  Five  header  lines  are  required  at  the 
beginning  of  the  data  entry.  These  lines  may  be  used  to 
describe  pertinent  information  about  the  loading  config¬ 
uration  to  be  analyzed.  This  information  will  be  piinted 
on  the  output  sheet  and  will  serve  to  identify  the  output. 
If  fewer  than  five  lines  are  used  to  identify  the  project, 
blank  lines  must  be  included  to  complete  the  required 

five  lines.  The  information  on  the  header  lines  may  be 
up  to  60  characters  maximum.  Data  for  example  1  for  the 
five  header  lines  are  as  follows: 

=SAMPLE  PROBLEM  FOR  SETTLEMENT 
=0CT  1978 

=VERTICAL  STRESS  DISTRIBUTION  INFLUENCE  FACTOR 
=UNIT  LOAD  OF  1.0  WITHOUT  FORCE  UNIT 
=TW0  RECTANGULAR  FOOTINGS 

b.  Loading  data.  The  type  and  number  of  lines  in  this  group 
vary  depending  upon  whether  stresses  are  from  footing 
loads,  an  embankment  load,  or  both  footing  and  embankment 
loads  are  being  analyzed.  The  type  of  loading  is  specified 
by  the  variable  RODE  described  below. 
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VOID  RATIO 


< 

^  0.700 
Q 


(0.850,  0.3) 


Cv  =  0.70  FT2/DAY 
K  =  0.0016  FT2/DAY 


Cc  =  CR  =  0.0543 


CR  =  0.0543 


(0.800,  2.5) 


Cv  =  0.70  FT2/DAY 
K  =  0.0016  FT2/DAY 


Cc  =  0.2491 


LOG  EFFECTIVE  STRESS  p  ,  KIPS/FT2 

Figure  5.  Void  ratio  versus  effective  stress  curve 
for  layer  2  in  Figure  3 


VOID  RATIO 


Table  2 

Pedometer  Test  Results  for  Clay  Layers  in  Example  Problem  1 


Layer 

No. 

Point 

No. 

Void 

Ratio 

Stress 
kips/f  t^ 

c 

c 

CR 

c 

V 

f  t^/day 

K 

f t^/day 

1 

1 

0.900 

0.2 

0.0500 

0.0500 

1.27 

0.0024 

1 

2 

0.850 

2 

0.3576 

0.0500 

1.27 

0.0024 

1 

3 

0.600 

10 

2 

1 

0.850 

0.3 

0.0543 

0.0543 

0.70 

0.0016 

2 

2 

0.800 

2.5 

0.2491 

0.0543 

0.70 

0.0016 

2 

3 

0.650 

10 

4 

1 

0.800 

0.2 

0.0402 

0.0402 

0.47 

0.0009 

4 

2 

0.750 

3.5 

0.3240 

0.0402 

0.47 

0.0009 

4 

3 

0.600 

10 

Table  3 

Results  of  Standard  Penetration  Test  of  Layer 
No.  3  (Sand  and  Gravel) 


Depth 

Corrected 

ft 

Blow  Count 

93.5  to  85 

32 

85  to  80 

35 

80  to  75 

20 

75  to  70 

30 

70  to  68.5 

40 

*  Average  blow  count  =  33.3. 
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(1)  The  first  line  in  this  group  is  the  following: 

KODE , NAREA 

(a)  Item  1 — KODE.  KODE  is  a  variable  which  indicates 
what  type  of  loading  configuration  is  to  be  used 
in  the  analysis.  If  KODE  is  input  as  1,  only 
uniform  rectangular  loads  are  to  be  used  in  the 
analysis.  If  KODE  is  input  as  2,  only  an  embank¬ 
ment  load  is  to  be  used.  If  KODE  is  entered  as  3, 
then  both  uniformly  loaded  rectangular  areas  and 
embankment  loads  are  to  be  used  in  the  analysis. 
(For  input  of  embankment  loads,  see  Appendix  A.) 

(b)  Item  2--NAREA.  NAREA  is  the  variable  indicating 
the  number  of  rectangular  uniformly  loaded  areas 
(footings)  to  be  entered.  NAREA  should  be  entered 
for  KODE  =  1  and  KODE  =  3  loading  conditions  but 
may  be  input  as  zero  for  KODE  =  2  (embankment 
only)  loading  conditions.  The  maximum  allowable 
value  of  NAREA  is  100. 

Input  for  example  1  for  this  data  line  is  as  follows: 

KODE, NAREA 

=1,  2 

(2)  The  next  line(s)  of  the  input  data  describes  the  loca¬ 
tion  and  loading  for  an  individual  rectangular  loaded 
area(s).  There  will  be  one  line  of  this  information 
for  each  rectangular  area  in  the  loading  configura- 
tion.  When  stresses  from  an  embankment  loading  only 
are  to  be  calculated  (KODE  =  2),  this  line  should  not 
be  included  in  the  input  data.  The  following  vari¬ 
ables  are  required  for  this: 

Q(I) , ZLAY (I),XC(1,I) ,YC(1,I) ,XC(2-4,I) , YC(2-4 , I) 

(a)  Item  1 — Q(I) ■  Q(I)  is  the  magnitude  of  the  uni¬ 
form  load  on  the  Ith  rectangular  area.  It  is  in 
units  of  LOAD/UNIT  AREA.  Any  units  for  weight 
and  length  may  be  used  as  long  as  all  input  data 
are  in  the  same  units. 

(b)  Item  2 — ZLAY(l).  Positive  ZLAY(I)  is  the  vertical 
distance  from  the  base  of  the  Ith  footing  to  the 
vertical  reference  plane  of  the  lowest  point  in 
the  embankment.  (No  footings  may  be  input  lower 
than  the  lowest  point  in  the  embankment.) 

(c)  Item  3 — XC(1 , I ) .  XC(1,I)  is  the  variable  name  of 

the  X  coordinate  of  the  first  corner  of  the  Ith 
rectangular  area.  The  dimensions  of  XC(1,I)  may 
be  in  any  units  compatible  with  the  remainder  of 
the  input  data. 
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(d)  Item  4 — YC(1, I) .  YC(1,I)  is  the  variable  name  of 

the  Y  coordinate  of  the  first  corner  of  the  Ith 
rectangular  area. 

(e)  Items  5-10— XC(2-4, 1)  and  YC(2-4,I).  These  are 
the  remaining  three  pairs  of  X  and  Y  coordinates 
which  define  the  corners  of  the  1th  rectangular 
area.  The  sides  of  the  area  do  not  have  to  be 
parallel  to  the  X  and  Y  axes,  but  the  corner 
points  should  be  input  in  either  clockwise  or 
counterclockwise  order  around  the  perimeter  of 
the  rectangular  area. 

Without  force  units  being  used  as  the  load,  this  would 
yield  a  factor  which  could  be  multiplied  times  any 
load  to  give  the  vertical  stress  at  that  point.  In¬ 
put  data  for  example  1  for  the  two  data  lines  (rectan¬ 
gular  areas  1  and  2,  respectively)  are  as  follows: 

Q(I) , ZLAY (I) ,XC(1,I) ,YC(1, I) ,XC(2, I) , YC(2, I) ,XC(3, I) , 
YC(3,I),XC(4,I),YC(4,I) 

=1.0  0.0  90.0  82.5  90.0  117.5  110.0  117.5  110.0  82.5 

Q(I)  ,  ZLAY  (1),XC(1,I),YC(1,I)  ,.\C(2 , 1 )  ,  Y(.  (2 , 1)  ,XC(J,I), 
YC(3,I),XC(4,I),YC(4,I) 

=1.0  0.0  120.0  120.0  120.0  135.0  130.0  135.0  130.0 

120.0 

£.  Stress  distribution.  This  group  of  data  defines  the  out¬ 
put  required  for  the  loading  conditions  described  in 
subparagraph  b  above.  The  output  may  be  in  two  forms, 
depending  on  the  needs  of  the  user.  The  first  type  of 
output  consists  of  values  of  vertical  stresses  printed 
along  a  vertical  line  in  the  X-Y-Z  plane.  For  this 
distribution,  the  values  of  X  and  Y  will  remain  constant. 
Stress  values  will  be  calculated  at  prescribed  increments 
between  prescribed  limits  along  the  vertical  line.  The 
second  type  of  output  option  consists  of  values  of  verti¬ 
cal  stresses  printed  at  increasing  values  of  X  along  a 
prescribed  line  in  the  X-Y  plane  at  a  constant  depth  (Z 
is  constant).  The  orientation  of  the  line  in  the  X-Y 
plane  is  defined  by  inputting  a  slope  and  an  intercept. 
There  is  no  limit  as  to  the  number  of  calculation  points 
on  a  given  distribution  or  on  how  many  distributions  may 
be  run  for  a  given  loading  configuration.  The  information 
for  a  single  stress  distribution  is  contained  on  two 
1 ines . 

(1)  The  input  variables  on  the  first  line  are: 

Xi!lST  ,  \W j-  M  ,  AMU 

(a)  Item  1 — NDIST.  NDIST  is  an  option  variable  which 
defines  whether  stress  distribution  in  a  vertical 
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or  horizontal  plane  is  required.  If  NDIST  is 
input  as  1,  a  vertical  plane  distribution  will  be 
assumed;  if  NDIST  is  input  as  2,  a  horizontal 
plane  distribution  is  calculated.  NDIST  also 
serves  to  indicate  when  all  the  stress  distribu¬ 
tions  for  a  given  loading  condition  are  completed. 
A  value  of  NDIST  equal  to  zero  will  cause  new 
header  cards  and  loading  data  to  be  read  in.  If 
no  new  loading  configuration  follows  NDIST  =  0, 
the  program  will  exit. 

(b)  Item  2 — NWEST.  NWEST  is  an  option  variable  which 
determines  whether  the  Westergaard  or  Boussinesq 
solution  will  be  used  to  determine  the  vertical 
stresses.  If  NWEST  =  0,  the  Westergaard  solution 
will  be  used;  if  NWEST  =  1,  the  Boussinesq  solu¬ 
tion  will  be  calculated. 

(c)  Item  3 — AMU.  AMU  represents  the  value  of 
Poisson's  ratio  to  be  used  in  the  Westergaard 
solution.  If  a  Boussinesq  solution  is  to  be  used 
(NWEST  =  I),  AMU  is  input  as  zero. 

Input  for  example  1  for  this  data  group  is  as  follows: 

NDIST,  NWEST,  AMT 

=2  1  0.0 

(2)  The  second  line  is  used  to  define  the  stress  distribu¬ 
tion  and  should  not  be  included  if  NDIST  =  0.  The 

card  includes  the  following  data: 

AINTL,  FINAL,  DELTA,  XP ,  YP ,  ZP ,  SLP  ,  BI.I NE 

Repeat  this  card  for  each  distribution  required 

(NDIST  times)  . 

(a)  Item  1 — AINTL.  AINTL  is  the  starting  point  coor¬ 
dinate  for  either  a  vertical  or  a  horizontal  plane 
distribution.  If  a  vertical  plane  distribution  is 
required  (NDIST  =  1),  the  value  of  AINTL  repre¬ 
sents  the  initial  (smallest)  depth  within  the 
range  of  the  distribution.  For  this  case,  AINTL 
must  be  positive.  In  the  case  of  a  horizontal 
plane  distribution  (NDIST  =  2),  AINTL  represents 
the  smallest  (initial)  X  coordinate  of  the  hori¬ 
zontal  pLane  distribution.  For  a  horizontal 
plane  distribution,  AINTL  may  be  positive  or 
negative. 

(b)  Item  2 — FINAL.  FINAL  is  t he  ending  point  coordi¬ 
nate  for  either  a  vertical  or  a  horizontal  plane 
distribution.  If  a  vertical  plane  distribution 

is  required  (NDIST  =  1),  the  value  of  FINAL  repre¬ 
sents  the  final  (largest)  depth  within  the  range 
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of  the  distribution.  For  this  case,  FINAL  must 
be  positive.  In  the  case  of  a  horizontal  plane 
distribution  (NDIST  =  2),  FINAL  represents  the 
largest  X  coordinate  of  the  horizontal  plane 
distribution.  For  this  case,  FINAL  may  be  posi¬ 
tive  or  negative. 

(c)  Item  3 — DELTA.  DELTA  is  the  distance  between 
calculation  points  for  both  a  horizontal  and  a 
vertical  plane  stress  distribution.  DELTA  should 
always  be  positive. 

(d)  Item  4 — XP.  XP  is  the  X  coordinate  for  the  loca¬ 
tion  of  stress  distribution  on  a  vertical  plane. 

If  stress  distribution  on  a  horizontal  plane  is 
required,  the  value  of  XP  may  be  input  as  zero. 

(e)  Item  5 — YP.  YP  is  the  Y  coordinate  for  the  loca¬ 
tion  of  a  vertical  plane.  If  stresses  on  a  hori¬ 
zontal  plane  are  required,  the  value  of  YP  may  be 
input  as  zero. 

(f)  Item  6 — ZP.  ZP  is  the  constant  depth  at  which  a 
horizontal  plane  is  located.  ZP  should  be  posi¬ 
tive  and  referenced  to  the  lowest  point  in  the 
embankment  or  footing  configuration.  If  a  verti¬ 
cal  plane  is  being  considered,  the  value  of  ZP 
may  be  input  as  zero. 

(g)  Item  7 — SLP.  SLP  is  the  slope  in  the  X-Y  plane 
of  the  line  along  which  stress  distribution  is 
required.  If  a  vertical  plane  is  being  consid¬ 
ered,  then  the  value  of  SLP  may  be  input  as  zero. 

(h)  Item  8 — BLINE.  BLINE  is  the  Y  intercept  of  the 
line  in  the  X-Y  plane  along  which  a  horizontal 
plane  stress  distribution  is  to  be  run.  The  value 
of  BLINE  may  be  input  as  zero  if  a  vertical  plane 
distribution  is  being  run. 

Input  for  example  1  for  this  data  group  is  as  follows: 

AINTL, FINAL, DELTA, XP , YP , ZP , SLP , BLINE 
=100.0  100.0  1.0  0.0  0.0  7.5  0.0  100.0 

21.  Input  data  for  example  1  as  stored  in  a  data  file  are  shown 
in  Table  4. 


22.  Output.  The  output  consists  of  vertical  stress  influence 
factors  at  the  points  requested.  For  example  1,  this  output,  along  with 
the  input  data,  is  shown  in  the  conversational  mode  in  Table  5.  Table  6 
shows  a  data  file  for  the  output  data. 
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Table  A 

Input  Data  File  for  Program  10016 
(Example  Problem  1) 


CPI9916  lSi  «t28  3/13/79 

liH  SAMPLE  PROBLEM  FOR  SETTLEMENT 
1019  OCT  1979 

1929  VERTICAL  STRESS  DISTRIBUTION  INFLUENCE  FACTOR 
1939  UNIT  LOAD  OF  1.9  UITHOUT  FORCE  UNIT 
1949  TWO  RECTANGULAR  FOOTINGS 
1959  1  2 

1969  1.9  9.9  99.9  82.5  99.9  117.5  119.9  117.5  119.9  82.5 
1979  1.9  9.9  129.9  129.9  129.9  13$. 9  139.9  13S.9  139.9  129.9 
1989  2  1  9.9 

1999  199.9  199.9  1.9  9.9  9.9  7.5  9.8  199.9 

1199  2  1  9.9 

1119  199.9  199.9  1.9  9.9  9.9  18.25  9.9  199.9 
1129  2  1  9.9 

1139  199.9  199.9  1.9  9.9  9.9  56.75  9.9  199.9 

1149  9  9  9.9 
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Table  5 


Input/Output  In  the  Conversational  Mode  for  Program  10016 
(Example  Problem  1) 


RUN  UESLIB/CORPS/I0016.R 

DO  VOU  WISH  TO  RUN  PROGRAM  FROM  EXISTING  DATA  FILE? 

•  h 

DO  VOU  UANT  OUTPUT  WRITTEN  TO  AN  OUTPUT  FILE? 

•  v 

INPUT  S  HEADER  LINES 
■  SAMPLE  PROBLEM  FOR  SETTLEMENT 
■OCT  1978 

•UERTICAL  STRESS  DISTRIBUTION  INFLUENCE  FACTOR 
•UNIT  LOAD  OF  1  0  UITHOUT  FORCE  UNIT 
•TUO  RECTANGULAR  FOOTINGS 
KODE.NAREA 
•1  2 


Om.ZLAY(I>.XCU.I>.YC(l,I).XC<2.n,YC(2.I).XCO.I),VC<3.I).XC(4.I),YC 

(4.1) 

•1  0  00  90  0  82  S  90  0  117  .5  110  0  117.5  110  0  82  5 

Q(I  ).ZLAY(I).XC(1.I  ).YC(1.I).XC<2.  I  ) . VC(2 . 1) , XC( 3 . 1 ) ,  YCO.  I ) ,  XC(  4  .  I ) ,  VC 
(4.1) 

■10  00  120  0  120  0  120  0  135  0  130  0  135  0  130  0  120  0 
NDIST .  NUEST,  AMU 
•2100 

AINTL  .  FINAL . DELTA. XP . VP , ZP . SLP . BLINE 
•100  0  100  0  10  00  00  75  0  0  100  0 

SAMPLE  PROBLEM  FOR  SETTLEMENT  PACKAGE 
OCT  1978 

UERTICAL  STRESS  DISTRIBUTION  INFLUENCE  FACTOR 
UNIT  LOAD  OF  1.1  UITHOUT  FORCE  UNITS 
TUO  RECTANGULAR  FOOTINGS 

BOUSSINESQ  SOLUTION 


HORIZONTAL  STRESS  DISTRIBUTION  AT  DEPTH(Z)  • 

ELASTIC  SOLUTION 

V-COORDINATE  X-COORDINATE  VERTICAL  STRESS 

100  00  100  09  9  882 


?  50 

NORMAL  LOADING 
UERTICAL  STRESS 

0  882 


NUMBER  OF  AREAS  USED  IN  CALCULATION  •  2 

NOTE-ALL  Z  VALUES  ARE  REFERENCED  TO  THE  LOWEST  PART  OF  THE  INPUT, 
CONFIGURATION 

NDIST.  NUEST.  AMU 
■2100 


AINTL .FINAL . DELTA . XP . VP , ZP . SLP .BLINE 
•109  0  100  0  1  0  0  0  0  0  18  25  0  0  100  0 


(Cont inued) 


Table  5  (Concluded) 


SAMPLE  PROBLEM  FOR  SETTLEMENT  PACKAGE 
OCT  1978 

VERTICAL  STRESS  DISTRIBUTION  INFLUENCE  FACTOR 
UNIT  LOAD  OF  18  UITHOUT  FORCE  UNITS 
TUO  RECTANGULAR  FOOTINGS 

BOUSSINESG  SOLUTION 

HORIZONTAL  STRESS  DISTRIBUTION  AT  DEPTH<2>  -  18. 25 

ELASTIC  SOLUTION  NORMAL  LOADING 
Y-COORDINATE  X-COORDINATE  VERTICAL  STRESS  VERTICAL  STRESS 


iee  ee  iee.ee  e.ses  e  se9 

NUMBER  OF  AREAS  USED  IN  CALCULATION  •  2 

NOTE-ALL  2  VALUES  ARE  REFERENCED  TO  THE  LOUEST  PART  OF  THE  INPUT. 
CONFIGURATION. 

NDIST .  NUEST.  AMU 

•2  i  e  e 

AINTL , FINAL . DELTA . XP , VP , ZP , SIP , BLINE 

•lee.e  lee.e  1  e  ee  e  e  S6  ?s  ee  iee  e 

SAMPLE  PROBLEM  FCR  SETTLEMENT  PACKAGE 
OCT  1978 

VERTICAL  STRESS  DISTRIBUTION  INFLUENCE  FACTOR 
UNIT  LOAD  OF  1.8  UITHOUT  FORCE  UNITS 
TUO  RECTANGULAR  FOOTINGS 

BOUSSINESG  SOLUTION 

HORIZONTAL  STRESS  DISTRIBUTION  AT  DEPTH(Z)  •  56  75 

ELASTIC  SOLUTION  NORMAL  LOADING 
V-COORDINATE  X-COORDINATE  VERTICAL  STRESS  VERTICAL  STRESS 


iee  ee  iee  ee  e  103  e  iea 

NUMBER  OF  AREAS  USED  IN  CALCULATION  >  2 

NOTE-ALL  Z  VALUES  ARE  REFERENCED  TO  THE  LOUEST  PART  OF  THE  INPUT. 
CONFIGURATION 

NDIST,  NUEST.  AMU 

•e  e  e.e 

* 
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Table  6 

Output  Data  File  for  Program  10016 
(Example  Problem  1) 


CCI8816  t4>59l  9  3/13/79 

SAMPLE  PRO* LEU  FOR  SETTLEMENT  PACKAGE 
OCT  1978 

VERTICAL  STRESS  DISTRIBUTION  INFLUENCE  FACTOR 
UNIT  LOAD  OF  1.8  WITHOUT  FORCE  UNIT 
TUO  RECTANGULAR  FOOTINGS 

BOUSSINESO  SOLUTION 

HORIZONTAL  STRESS  DISTRIBUTION  AT  DEPTH(Z)  •  7, SB 

ELASTIC  SOLUTION  NORnAL  LOADING 
V-COORDINATE  X-COORDINATE  VERTICAL  STRESS  VERTICAL  STRESS 


1BB.BB  188.88  8.882  8.882 

NUMBER  OF  AREAS  USED  IN  CALCULATION  •  2 

NOTE-ALL  Z  VALUES  ARE  REFERENCED  TO  THE  LOWEST  PART  OF  THE  INPUT, 
CONFIGURATION. 

SAMPLE  PROBLEM  FOR  SETTLEMENT  PACKAGE 
OCT  1978 

VERTICAL  STRESS  DISTRIBUTION  INFLUENCE  FACTOR 
UNIT  LOAD  OF  1.8  WITHOUT  FORCE  UNIT 
TWO  RECTANGULAR  FOOTINGS 

BOUSSINESO  SOLUTION 

HORIZONTAL  STRESS  DISTRIBUTION  AT  DEPTH (Z>  •  18.25 

ELASTIC  SOLUTION  NORNAL  LOADING 
V-COORDINATE  X-COORDINATE  VERTICAL  STRESS  VERTICAL  STRESS 


*••.80  188.88  8.589  8.589 


NUMBER  OF  AREAS  USED  IN  CALCULATION  •  2 

NOTE-ALL  Z  VALUES  ARE  REFERENCED  TO  THE  LOWEST  PART  OF  THE  INPUT, 
CONFIGURATION. 

SAMPLE  PROBLEM  FOR  SETTLEMENT  PACKAGE 
OCT  1978 

VERTICAL  STRESS  DISTRIBUTION  INFLUENCE  FACTOR 
UNIT  LOAD  OF  1 .8  WITHOUT  FORCE  UNIT 
TUO  RECTANGULAR  FOOTINGS 

BOUSSINESO  SOLUTION 

HORIZONTAL  STRESS  DISTRIBUTION  AT  DEPTH(Z)  •  SG.7S 

ELASTIC  SOLUTION  NORMAL  LOADING 
V-COORDINATE  X-COORDINATE  VERTICAL  STRESS  VERTICAL  STRESS 


188-88  188.88  8.183  8.183 

NUMBER  OF  AREAS  USED  IN  CALCULATION  •  2 

NOTE -ALL  Z  VALUES  ARE  REFERENCED  TO  THE  LOWEST  PART  OF  THE  INPUT, 
CONFIGURATION. 
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Program  MAGS ET I I 

23.  Organization  of  input.  Input  for  MAGSETII  is  broken  down  into 
two  basic  areas:  problem  control  and  data  entry.  The  first  governs  the 
execution  of  the  program  and  allows  the  user  to  describe  what  data  and 
what  form  of  data  are  to  be  entered.  One  more  option  to  the  program  has 
been  added  to  perform  a  rate  of  settlement  analysis. 

24.  Mode  of  input.  Input  to  the  program  can  be  either  from  the 

terminal  or  from  a  data  file.  All  input  is  in  free  field.  Data  items 

can  be  separated  by  a  blank  or  a  comma.  If  the  program  is  being  run 

from  a  data  file,  lines  of  the  data  file  must  be  preceded  by  a  line 
number.  When  the  program  is  run  from  a  terminal,  it  can  create  files  to 
save  the  input  data  and  output  from  the  run. 

25.  After  receiving  the  output  from  10016,  MAGSETII  can  be  used 

to  calculate  the  settlement  and  the  rate  of  consolidation.  Input  for 

this  program  is  shown  in  Tables  7  and  8. 

26.  Problem  control  input.  The  first  line  in  this  section  gives 
the  information  on  one  particular  run  and  controls  whether  or  not  more 
than  one  problem  is  going  to  be  run.  The  next  line  is  a  title  descrip¬ 
tion  of  the  particular  problem. 

27.  Data  input.  The  first  line  in  this  section  gives  the  problem 
options  for  the  output  and  data  input  control.  There  are  eight  of  these 
options : 

<i.  Unit  indicator.  This  specifies  whether  the  units  are  to 
be  shown  in  the  output. 

_b.  Effective  stress  indicator.  This  indicates  whether  the 
effective  stress  is  to  be  calculated  at  midpoints  or 
input  at  each  soil  layer  or  to  combine  the  calculated  and 
the  input  effective  stress. 

£.  Effective  stress  history  specifications.  This  indicates, 
if  clay  layers  are  present,  whether  the  effective  stress 
history  is  to  be  input  for  each  clay  layer  or  one  is  to 
be  used  for  all  clay  layers  or  the  vertical  stress  distri¬ 
bution  function  is  to  be  multiplied  time  one  effective 
stress  history  for  all  clay  layers. 

d.  Deformation  curve  type.  This  indicates,  if  clay  layers 
are  present,  whether  the  deformation  curve  is  a  strain 
or  a  void  ratio  versus  effective  stress  relationship. 
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Table  7 

Input  in  the  Conversational  Mode  for  Program  MAGSETII 
(Example  Problem  1) 


RUN 


INPUT  NAME  OP  DATA  FILE.  HIT  A  CARRIAQE  RETURN 
If  DATA  IS  TO  SE  READ  FRON  THE  TERNINAL. 


INPUT  A  FILE  NAME  FOR  DATA  IN  g  CHARACTERS  OR  LESS 

HIT  A  CARRIAOE  RETURN  IF  YOU  DO  NO  UANT  TO  SAVE  THIS  FILE 

•REEDIN 

INPUT  A  FILE  NAME  FOR  OUTPUT  IN  8  CHARACTERS  OR  LESS 

HIT  A  CARRIAQE  RETURN  IF  DATA  IS  TO  BE  PRINTED  ON  TERMINAL 

-REEDOUT 


INPUT  PROBLEM  CONTROL  INFORMATION 
NPROB  •  PROBLEM  NUMBER 

NLAVER  -  NUMBER  OF  SOIL  LAVERS  IN  PROFILE  HAMMS 
NLAST  -  •  IF  CURRENT  PROBLEM  ISN'T  LAST  ONE 

1  CURRENT  PROBLEM  IS  LAST  ONE  IN  DATA  SET 

M.S.l 


INPUT  PROBLEM  OPTIONS 

I OPT  < l )  -  UNITS  INDICATOR 

1  -  UNITS  TO  BE  PRINTED  IN  OUTPUT 

2  -  UNITS  KILL  NOT  BE  PRINTED 
IOPT(g)  -  INSITU  EFFECTIVE  STRESS  INDICATOR 

1  -  IES  CALCULATED  AT  MIDPOINTS 
8  -  IES  INPUT  AT  EACH  LAVER 

3  -  IES  INPUT  AT  EACH  LAVER  AND  ADDED  TO  CALC. IES 
I  OPTO)  -  EFFECTIVE  STRESS  HISTORY  SPECIFICATIONS 

S  -  NO  CLAV  LAVERS 

1  -  ESH  INPUT  FOR  EACH  CLAV  LAVER 

2  -  ONE  ESH  INPUT  AND  USED  FOR  ALL  LAVERS 

3  -  ONE  ESH  INPUT  AND  USED  FOR  ALL  CLAV  LAVERS 

STRESS  DISTRIBUTION  FUNCTION  UILL  BE  INPUT 
I  OPTO )  -  DEFORMATION  CURVE  TYPE 
B  -  NO  CLAV  LAVERS 

1  -  STRAIN-EFFECTIVE  STRESS  CURVES 

2  -  VOID  RATIO-EFFECTIVE  STRESS  CURVES 
IOPT(S)  -  DEFORMATION  CURVE  SPECIFICATION 

•  -  NO  CLAV  LAVERS 

1  -  DC  INPUT  USINQ  COORD.  PTS.  OF  VOID  RAIO  OR 

STRAIN  VS.  EFFECTIVE  STRESS 

2  -  DC  INPUT  USINQ  SLOPES. EFF  STRESS  VALUES  AND  A 

REFERENCE  COORDINATE 

IOPT(S)  -  SAND  SETTLEMENT  METHOD  INDICATOR 
B  -  NO  SAND  LAVERS 
1  •  HEVERHOFF'S  METHOD 
8  -  D'APPOLONIA'S  METHOD 

3  -  ALL  THREE  METHODS 

IOPT<?)  -  VERTICAL  STRAIN  INFLUENCE  FUNCION 
B  -  NO  SAND  LAVERS 
1  -  CURVE  Q 
8  -  CURVE  F 

3  -  VERTICAL  STRAIN  INFLUENCE  FUNCION  UILL  BE  INPUT 
IOPT(S)  -  DATUM  CONVERSION  OPTION 

(Continued) 
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Table  7  (Continued) 


INPUT  TITLE 

TITLE  -  DESCRIPTION  OF  PROILER  IN  SS  CHARACTERS  OR  LESS 
•  SAMPLE  PROILEN  FOR  SETTLEMENT  PACKAGE 


1  -  DEPTHS 
LL  NOT  »E  CONCERTED 

R  -  DEPTHS 
LL  IE  CONVERTED 
•1.1. 3. 8. 1,1, 1.8 


OR  ELEVATIONS  HI 
OR  ELEVATIONS  UI 


INPUT  UNITS 

lUNIT(l)  -  LENQTH  UNITS  IN  COLUMNS  1-1« 
I  UNIT ( 8 )  -  FORCE  UNITS  IN  COLUMNS  17-38 
•  FEET  KIPS 


INPUT  GROUND  HATER  DATA 

OH  -  TNIT  HEIGHT  OF  HATER 

QHELEV  -  DEPTH  OR  ELEU  OF  GROUND  HATER  SURFACE 
•  G684.1GG  G 


INPUT  LAVER  INTERFACE  INFORMATION 
L  -  LAVER  INTERFACE  NUNIER 
DEPTH! L)  -  DEPTH  OR  ElEV  TO  TOP  OF  LAVER 

LAVER  NUMBER  1 
•1.185  G 

LAVER  NUMBER  8 
•8. 115  G 

LAVER  NUMBER  3 
•3/1GG.G 

LAVER  NUMBER  4 
•4,03. 5 

LAVER  NUMBER  6 
■S.SS.S 

LAVER  HUMBER  S 
•6.48  G 


(Continued) 
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Table  7  (Continued) 


INPUT  SOU  PROPERTIES 
L  -  LAVER  NUMBER 
NTVPf <L>  -  SOIL  TYPE 
1  -  CLAV 
t  -  SANS 

3  -  INCOMPRESSIBLE 
WOTIl)  -  TOTAL  UNIT  WEIGHT  OP  SOIL 

LAVER  NUMBER  1 
•1.3. .ISO 

LAYER  NUMBER  S 

•8.1.. 1«S 

LAVER  NUMBER  3 
•3.1.  ISO 

LAVER  NUMBER  4 
•4.8. .130 

LAVER  NUMBER  S 

•5.1.. 130 

INPUT  DATUM  CONVERSION  INFORMATION 
DATUM  -  DATUM  ELEVATION 

DIFELV  -  DIFF.  IN  ELEV.  BETWEEN  DATUM  ELEV  1  TOP  OF  1ST  LAVER 
•185.S.S.S 

INPUT  EFFECTIVE  STRESS  INCREMENTS 
SIQKNS)  -  THE  < NSTH )  ESI 
LS  -  LAST  INCREMENT  INDICATOR 

0  -  IF  MORE  ESI'S  TO  BE  INPUT 
1  -  IF  LAST  ESI 

STRESS  INCREMENT  NUMBER  1 

— 1  0,0 

STRESS  INCREMENT  NUMBER  8 
•1  3.0 

STRESS  INCREMENT  NUMBER  3 

•a  o.o 

STRESS  INCREMENT  NUMBER  4 

•a  s.o 

STRESS  INCREMENT  NUMBER  5 
— I  5.1 

INPUT  STRESS  DISTRIBUTION  FUNCTION 
L  -  CLAV  LAVER  NUMBER 
F(L>  -  VALUE  OF  SDF 

LAVER  NUMBER  Z 

•a.  «« 

LAVER  NUMBER  3 
•3. .51 

LAVER  NUMBER  5 
•5.  103 


(Continued) 
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Table  7  (Continued) 


INPUT  DEFORMATION  CURVE  COORDINATE  POINTS 
I LAVER  -  CLAV  LAVER  NUNIER 
L INEPT ( 1 , IPT)  -  FIRST  POINT  ON  DC  •  l 
Ed.  IPT)  -  VOID  RA10  OR  STRAIN  COORD  AT  1ST  PT  ON  DC 
SIONAPd.IPT)  -  EFFECTIVE  STRESS  COORD  AT  1ST  PT  ON  DC 
•2.1  .  .  9S.  81 


INPUT  DEFORMATION  CURVE  -  SUSSEOUENT  COORD  PTS 
I LAVER  -  CLAV  LAVER  NUNIER 
LINEPT(I.IPT)  -  COORD  PT  ON  DC 
Ed. IPT)  -  VOID  RATIO  OR  STRAIN  AT  COORD  PT 
SIONAPd.IPT)  -  EFFECTIVE  STRESS  AT  COORD  PT 
ERd. IPT-1)  -  VOID  RATIO  OR  STRAIN  COORD  TO  SE  USED  TO 
CALCULATE  EXPANSION  SLOPE 

SIQRd.  IPT-1 )  -  EFFECTIVE  STRESS  TO  IE  USED  TO  CALC 
EXPANSION  SLOPE 
IP  -  LAST  POINT  INDICATOR 
S  -  NOT  LAST  POINT 
1  '  LAST  POINT 
INCREMENT  I 
•2.2.  15.2  .0.  90. .20.0 
INCREMENT  9 
•2.3.  60.10  0. .85. 2. 0.1 

INPUT  DEFORMATION  CURVE  COORDINATE  POINTS 
UAVER  -  CLAV  LAVER  NUNIER 
L INEPT d . IPT )  -  FIRST  POINT  ON  DC  •  1 
E(I. IPT)  -  VOID  RAIO  OR  STRAIN  COORD  AT  1ST  PT  ON  DC 
SIONAPd.IPT)  -  EFFECTIVE  STRESS  COORD  AT  1ST  PT  ON  DC 
•3.1.  If.. 30 


INPUT  DEFORMATION  CURVE  -  SUISEOUENT  COORD  PTS 
UAVER  -  CLAV  LAVER  NUNIER 
LlNEPTd.IPT)  -  COORD  PT  ON  DC 
Ed. IPT)  -  UOID  RATIO  OR  STRAIN  AT  COORD  PT 
SIONAPd.IPT)  -  EFFECTIVE  STRESS  AT  COORD  PT 
ERd.  IPT-1)  -  VOID  RATIO  OR  STRAIN  COORD  TO  IE  USED  TO 
CALCULATE  EXPANSION  SLOPE 

SIORd.IPT-l)  -  EFFECTIVE  STRESS  TO  IE  USED  TO  CALC 
EXPANSION  SLOPE 
LP  -  LAST  POINT  INDICATOR 
0  -  NOT  LAST  POINT 
1  -  LAST  POINT 
INCREMENT  8 
•3.2.  IS. 2  56. .16,  30.0 
INCREMENT  3 
•3.3.  65.10  0.  10.8  6.1 


(Continued) 


Table  7  (Continued) 


INPUT  DEFORMATION  CURVE  COORDINATE  POINT* 

I LAVER  -  CLAV  LAVER  NUMBER 
UNEPTt I. IPT )  -  FIRBT  POINT  ON  DC  ■  1 
EC  I. IPT)  -  VOID  RAIO  OR  STRAIN  COORD  AT  1ST  PT  ON  DC 
SIQNAPC I . IPT  >  -  EFFECTIVE  STRESS  COORD  AT  1ST  PT  ON  DC 
•5,1.  SO.  BO 

INPUT  DEFORMATION  CURVE  -  SUBSEQUENT  COORD  PTS 
ILAVER  -  CLAV  LAVER  NUMBER 
UNEPTt  I.  IPT)  -  COORD  PT  ON  DC 
EU.1PT)  -  VOID  RATIO  OR  STRAIN  AT  COORD  PT 
SIQNAPt I . IPT )  -  EFFECTIVE  STRESS  AT  COORD  PT 
ERU.IPT-l)  -  VOID  RATIO  OR  STRAIN  COORD  TO  BE  USED  TO 
CALCULATE  EXPANSION  SLOPE 

SICRCI.IPT-l)  -  EFFECTIVE  STRESS  TO  BE  USED  TO  CALC 
EXPANSION  SLOPE 
LP  -  LAST  POINT  INDICATOR 
S  -  NOT  LAST  POINT 
1  -  LAST  POINT 
INCREMENT  2 
•S.2.  75,3.51,  SO.  20, S 
INCRENENT  3 

•S.3.  60.10  0. .75.3  SO.l 


INPUT  PENETRATION  RESISTANCE 
L  -  SAND  LAVER  NUMBER 

BLOU(L)  -  STANDARD  PENETRATION  TEST  BLOUCOUNT  IN 
BLOUS  PER  FOOT 
0  -  IF  I0PT(S)-3 

OCtL)  -  STATIC  CONE  PENETRATION  RESISTANCE 
0  -  IF  I0PT(6).l  OR  2 

SAND  LAVER  NUMBER  4 
•4.33  3.0  0 


INPUT  FOOTINO  DATA 

FP  -  AUG  FOOTING  PRESSURE 
FB  -  FOOTING  UIDTH 
F DEPTH  -  DEPTH  OR  ELEV  OF  FOOTINO 
•5  0.20  0.115  0 

INPUT  NEVERHOFF'S  CONVERSION  FACTORS 
CONUFT  -  LENGTH  CONVERSION  FACTOR 
CONUTN  -  FORCE  CONVERSION  FACTOR 
•1  0.0  5 


(Continued) 


Table  7  (Concluded) 


INPUT  PROGRAM  CONTROL 

NHIST  -  0  NO  RATE  OF  CONSOLIDATION  IN  OUTPUT 

1  HISTORV  OF  RATE  OF  CONSOLIDATION  IN  OUTPUT 
8  SETTLEMENTS  AND  DEGREE  OF  CONSOLIDATION  AT 
THE  END  OF  LOADING  INCREMENTS  AND  SPECIFIC 
TIMES  AFTER  LOADIND 

•2 


INPUT  COEFFICIENT  OF  CONSOLIDATION 
L  -  LAVER  NUMBER 

CU  -  COEFFICIENT  OF  CONSOLIDATION! SO  FT. /DAY) 

NTOP  -  S  IF  TOP  IS  FREELV  DRAINED 

1  IF  TOP  IS  NOT  FREELV  DRAINED 
NDOT  -  *  IF  BOTTOM  IS  FREELV  DRAINED 

1  IF  BOTTOM  IS  NOT  FREELV  DRAINED 

CU  FOR  LAVER  NUMBER  2 
■2.1  27. S. I 

CU  FOR  LAVER  NUMBER  3 
•3.  70.1.# 

CU  FOR  LAVER  NUMBER  S 
•5. .47.8,8 

INPUT  TIMES  FOR  LOADIND  INCREMENTS 

T  -  THE  TIME! IN  DAYS)  FROM  THE  BfQINNINQ  OF 

CONSTRUCTION  TO  THE  END  OF  THE  LOAD  INCREMENT 
!LOAD  INCREMENT ‘STRESS  INCREMENT  IN  OPTION  3) 

MAX.  NUMBER  •  18 

LOAD  INCREMENT  NO  1 

•50  0 

LOAD  INCREMENT  NO.  2 
-75  8 

LOAD  INCREMENT  NO.  3 

•288  8 

LOAD  INCREMENT  NO  4 

•388  8 

LOAD  INCREMENT  NO  5 
•358  8 

INPUT  TIMES  FOR  AFTER  CONSTRUCTION 

T  -  NUMBER  OF  DAYS  FROM  THE  BEOINNINQ  OF  CONSTRUCTION 
TO  A  TIME  AFTER  THE  FINAL  LOADIND 
ENTER  0.0  FOR  LAST  ONE. ! MAX  NUMBER* 1#) 

A  TIME  AFTER  CONSTRUCTION  1 

•400.0 

A  TIME  AFTER  CONSTRUCTION  • 

•588  8 

A  TIME  AFTER  CONSTRUCTION  3 

•688.8 

A  TIME  AFTER  CONSTRUCTION  4 

•1888  8 

A  TINE  AFTER  CONSTRUCTION  6 

*0  0 


t 
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Table  8 

Input  Data  File  for  Program  MAGS KT 11 
(Example  Problem  1) 


REEDIN  8  43  t4  3/ 


10600 

1 

5 

1 

10010 

STAPLE 

PROBLEM  FOR 

SETTLEMENT 

PACKAGE 

10020 

1 

1 

3  2 

1  1 

1  2 

10030 

FEET 

KIPS 

10040 

0 

0624 

100  0000 

10050 

1 

125 

0000 

10060 

2 

115 

0000 

10070 

3 

100 

0000 

10080 

4 

93 

5000 

10090 

5 

68 

5000 

10100 

6 

48 

0000 

10110 

1 

3 

0.1000 

10120 

2 

1 

0  1200 

10130 

3 

1 

0  1200 

10140 

4 

2 

0  1300 

10150 

5 

1 

0  1300 

10160 

125 

0000 

0 

10170 

-1 . 

0000 

0 

10180 

1 

3000 

0 

10190 

2 

0000 

0 

10200 

2 

5000 

0 

10210 

-0 

5000 

l 

10220 

2 

0 

8800 

10230 

3 

0 

5100 

10240 

5 

0 

1030 

10250 

2 

1 

0  9000 

0  2000 

10260 

2 

2 

0  8500 

2  0000 

0  9000 

0  2000 

0 

10270 

2 

3 

0  6000 

10  0000 

0  8500 

2  0000 

1 

10280 

3 

1 

0  8500 

0  3000 

10290 

3 

2 

0  8000 

2  5000 

0  8500 

0  3000 

0 

10300 

3 

3 

0  6500 

10  0000 

0  8000 

2  5000 

1 

10310 

5 

1 

0  8000 

0  2000 

10320 

5 

2 

0  7500 

3  5000 

0  8000 

0  2000 

0 

10330 

5 

3 

0  6000 

10  0000 

0  ?500 

3  5000 

1 

10340 

4 

33 

3000  0. 

10350 

5 

0000 

200006 

115  0000 

10360 

1 

0000 

0  5000 

10380 

2 

10390 

2 

1 

2700  0 

1 

10400 

3 

i 

1  7000  1 

0 

10410 

5 

« 

i  4700  0 

0 

10420 

50  00 

10421 

75  00 

18422  288  80 

10423  300  80 

18424  350  00 

10425  400  00 

10426  500  00 

1042?  600  00 

10428  1000  0 

10429  0 
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e.  Deformation  curve  specifications.  This  indicates  if  the 
deformation  curves  are  to  be  input  by  coordinate  points 
upon  a  void  ratio  or  strain  versus  effective  stress  curve 
with  slopes  co  be  calculated  between  points  entered  or  by 
entering  the  slopes  and  reference  points  on  the  curve. 

jf.  Sand  settlement  methods  indicator.  If  a  sand  layer  is 
present,  it  indicates  the  method  of  analysis  that  needs 
to  be  employed  by  the  program  for  estimating  the  settle¬ 
ment  of  the  sand  layers. 

£.  Vertical  strain  influence  functions.  If  sand  layers  are 
present,  it  indicates  whether  one  of  the  built-in  strain 
influence  curves  is  to  be  used  or  a  function  is  to  be 
entered  by  the  user. 

_h,  Datum  conversion  option.  This  is  used  to  select  whether 
the  depth  or  elevation  is  to  be  converted  for  the  output. 

28.  Output .  The  output  from  a  successfully  executed  MAGSETII 
problem  is  printed  under  several  headings.  The  information  under  these 
headings  may  vary  slightly,  depending  on  the  problem  options  chosen.  A 
brief  description  of  the  information  printed  under  these  headings  is 
given  below. 

a.  Problem  specifications.  Printed  under  this  heading  are 
the  program  header,  title,  and  units. 

ja.  Soil  profile  description.  The  soil  profile  description 
prints  the  layer  number,  layer  type,  interface  depths  or 
elevations,  datum  elevations,  layer  thickness,  cad  the 
total  unit  weights  of  the  soils.  Also  under  this  heading 
are  the  groundwater  information,  the  unit  weight  of  water, 
the  depth  or  elevation  of  the  groundwater  table,  and  the 
datum  elevation  of  the  groundwater. 

cL.  In  situ  effective  stress.  Under  this  heading  are  the 
input  and  in  situ  effective  stresses  in  each  layer  and 
the  in  situ  effective  stress  used  by  the  program. 

d.  Clay  settlement  data.  This  section  contains  the  effec¬ 
tive  stress  history  and  deformation  curve  data.  The  in¬ 
put  effective  stress  increments  and  the  effective  stress 
history  are  printed.  The  input  data  used  to  specify  the 
deformation  curves  are  printed  along  with  any  data  calcu¬ 
lated  that  define  the  deformation  curve.  If  void  ratio 
versus  effective  stress  curves  are  input,  the  compression 
and  expansion  indexes  C  and  C  are  output  along  with 
with  strain  compression  and  strain  expansion  indexes  C 
and  C  .  If  strain  versus  effective  stress  curves  are 
input,  only  the  strain  compression  and  strain  expansion 
slopes  are  output. 


35 


ia.  Sand  settlement:  data.  This  section  contains  the  data 
used  in  the  sand  settlement  calculations.  It  includes 
calculation  methods,  foundation  data,  the  penetration 
resistance  for  each  sand  layer,  and  the  strain  influence 
function  used.  It  also  includes  information  which  is 
method-dependent,  such  as  D' Appolonia' s  parameters  and 
Meyerhof's  conversion  factors. 

_f.  Clay  settlement  contributions.  The  clay  settlement  con¬ 
tributions  contain  the  settlement  in  each  layer  due  to 
each  effective  stress  increment,  the  total  clay  settle¬ 
ment  in  each  layer,  the  settlement  in  the  clay  profile 
due  to  each  effective  stress  increment,  and  the  total 
clay  settlement. 

j>.  Clay  compressibilities.  The  coefficient  of  constrained 
compressibility  m  in  each  layer  for  each  effective 
stress  increment  is  printed.  The  column  header  El  defines 
the  void  ratio  or  strain  value  at  the  beginning  of  the 
effective  stress  increment  depending  on  the  form  of  the 
deformation  curves  input.  The  column  header  E2  defines 
the  void  ratio  or  strain  at  the  end  of  the  effective 
stress  increment.  The  column  header  DELTA  E  i  •  the  value 
of  E2  minus  El. 

h.  Sand  settlements.  The  settlements  in  each  sand  layer 
are  printed  under  the  method  of  analysis.  The  total 
sand  settlement  over  the  sand  profile  is  printed  along 
with  the  total  clay  settlement  and  total  profile 
settlement. 

iL.  Error  messages.  Various  checks  are  made  on  the  input 
data.  If  any  of  these  checks  fail,  an  error  message  is 
printed  and  the  program  terminates  execution.  The  error 
messages  have  been  worded  to  be  reasonably  self- 
explanatory.  If  confusion  results,  however,  refer  to 
Chapter  1  of  Schiffman,  duber.ville,  and  Partyka  (1976) 
under  the  appropriate  sections.  Table  9  presents  the 
output  of  MAGSETII  for  example  problem  1. 

j_.  Degree  of  consolidation.  The  time  (in  days),  time  factor 
(TV),  and  the  degree  of  consolidation  are  output  for  each 
soil  layer  at  each  10  percent  increment  or  at  a  specific 
time. 

Comparison  with  hand  calculations 

29.  Example  problem  1  was  worked  by  hand  using  conventional 
methods.  Results  from  program  10016,  using  a  Boussinesq  solution,  were 


compared  to  answers  from  an  influence  chart  for  vertical  pressure  for 
Boussinesq's  equation,  commonly  known  as  Newmark's  chart 


Table  9 

Output  Data  for  Program  MAGSKTII 
(Example  Problem  1) 


REEDOUT 


1 5 1 26 1 1 S  2^13x79 


t 

t  NAGSET-II 

t 

t  MAGNITUDE  OF  SETTLEMENT  OF 
t  A  HULTI -LAYERED  SOU  SYSTEM 
t 

iimitiitKitififitftfttttm 


ttttttttttttttttttttst 

t  SPECIFICATIONS  FORt 
S  PROBLEM  NO.  1  » 

imtutmmmuut 

tint  TITLE  ttstt 

SAMPLE  PR01LEH  FOR  SETTLEMENT  PACKAGE 

ttttt  UNITS  ttttt 
LENGTH  FORCE 

FEET  KIPS 

tttttttltlltllll 

*  SOIL  PROFILE  * 
t  DESCRIPTION  t 

ttttiititmmi 


DATUM  ELEVATION 

125.00 

DIFFERENCE 

IN  ELEVATION  • 

B. 

LAYER 

SOIL 

INTERFACE 

DATUM 

UNIT 

NUMBER 

TYPE 

ELEVATION 

ELEVATIONS 

THICKNESS 

WEIGHT 

12S.M 

125. 06 

1 

INCOMP 

115. M 

115. •• 

1B.BB 

•  .IBM 

2 

CLAY 

HIM 

1M.M 

1S.M 

B.12M 

3 

CLAY 

93.  SB 

93.50 

6. SB 

B.12M 

4 

SAND 

68. SB 

68. SB 

25.  BO 

0.13M 

S 

CLAY 

48. M 

48. M 

2B.SB 

B.13BB 

UNIT  WEIGHT 

GROUND  WATER  GROUND 

WATER 

OF  WATER 

LEVEL 

DATUM  ELEVATION 

•  •624 

1M.BB 

1M. 

•• 

iimttmttititmmmti 
S  INSITU  EFFECTIVE  STRESS  * 
ittittttttttftiiiiimtmt 


LAYER 

INPUT 

CALCULATED 

NUMBER 

VALUE 

VALUE 

1 

- 

B.SBBB 

2 

- 

1.9BBB 

3 

- 

2.9872 

4 

- 

4.B194 

5 

- 

5. 5573 

INSITU 

STRESS 

•  .SOM 

1.9SH 

2.9872 

4.8194 

S.5S73 


ttmttttitmmtttittt 
t  CLAV  SETTLEMENT  DATA  t 
ttttutttitttmttstiiti 


EFFECTIVE  STRESS  INCREMENTS  INPUT  IV 

A  DISTRIBUTION  FUNCTION 
ttttt  STRESS  INCREMENTS  FOR  STRATUM  ttttt 
POINT  NUMBER  STRESS  INCREMENT 


-1 


•••• 

1.3IM 

2.0*09 

2. SMB 

-B.SMB 


(Continued) 
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Table  9  (Continued) 


REEDOUT 

1S>26> 15 

2/13/79 

ttttt  STRESS 

DISTRIBUTION  FUNCTION  ttttt 

LAYER 

NUMBER 

VALUE 

2 

M8M 

3 

•.site 

S 

•.1*36 

ttttt  EFFECTIVE  STRESS  HISTORY  ttttt 

LAVER  HO. 

PT.  HO. 

STRESS  INCREMENT  STRESS  VALUES 

2 

1 

-0 . 8890 

1.0200 

2 

2 

1.1440 

2.1640 

2 

3 

1 . 7600 

3.9240 

2 

4 

2.2000 

6.1240 

2 

S 

-0.4400 

5.6840 

3 

1 

-0.5100 

2.4772 

3 

2 

0.6630 

3.1402 

3 

3 

1.0200 

4.1602 

3 

4 

1 . 2750 

5.4352 

3 

S 

-0.2550 

5.1802 

5 

1 

-0.1030 

S.4543 

5 

2 

0. 1339 

5.5882 

5 

3 

0.2060 

S.7942 

5 

A 

0.257S 

6.0517 

5 

S 

-0.0515 

6.0002 

DEFORMATION  CURVES  INPUT  RV 

♦ 

COORDINATE  POINTS 

ttttt  COORDINATES  Of 

POINTS  ON  THE  DEFORMATION 

CURVES  ttttt 

LAYER 

POINT 

REBOUND 

REBOUND 

NUMBER 

NUMBER 

VOID  RATIO 

STRESS 

VOID  RATZO 

STRESS 

2 

1 

0.9000 

0.2000 

2 

2 

O.BSOO 

2.0000 

0.9000 

0.2000 

2 

3 

0.6000 

10.0000 

0.1500 

2.0000 

3 

1 

0 . 8500 

0.3000 

3 

2 

0.8000 

2.5000 

0.8500 

0.3000 

3 

3 

0.6500 

10.0000 

0.8000 

2.5000 

5 

1 

0.8000 

0.2000 

5 

2 

0.7500 

3.S000 

0.8000 

0.2000 

5 

3 

0.6000 

10.0000 

0.7500 

3.5000 

ttttt  SLOPES 

ON  THE  DEFORMATION  CURVES  ttttt 

LAVER 

LINE 

CC 

CE 

CC 

CE 

NUMBER 

NUMBER 

(STRAIN) 

(STRAIN) 

2 

1 

0.0500 

0.0500 

0.0263 

0.0263 

2 

2 

0.3577 

0.3577 

0.1933 

0.1933 

3 

1 

0.0543 

0.0543 

0.0294 

0.0294 

3 

2 

0.2491 

0.2491 

0.1384 

0.1384 

6 

1 

0.0402 

0.0402 

0.0223 

0.0223 

5 

2 

0.3290 

0.3290 

0.1880 

0.1880 

(Continued) 
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Table  9  (Continued) 


*  SAND  SETTLEMENT  DATA  t 


(MM  FOUNDATION  DATA  ***** 


SUiMHJT 

15*261 IS  2/13/79 

FOOTING 

FOOTING 

FOOTING 

uipth 

PftCSSUftC 

♦ 

ELEVATION 

20.  MOO 

5.0000 

115.0000 

Util  SAND  SETTLEMENT  METHODS  mil 
NEVERHOF 


*****  PENETRATION  RESISTANCE  1***1 
LAVER  SPT  STATIC  CONE 

NUMDER  1LOUCOUNT  PENETRATION  RESISTANCE 

4  33.3* 


*****  MEYERHOF  UNITS  CONVERSION  FACTORS  ***** 

one  length  unit  eouals  i.rro*  feet 

ONE  FORCE  UNIT  EOUALS  R.SRRR  TONS 


ttttt  STRAIN  INFLUENCE  FUNCTION  ***** 

THE  STRAIN  INFLUENCE  FUNCTION  USED  IS  CURVE  C 


I  CLAV  SETTLEMENT  COMTH  HUT  IONS  t 


*****  SETTLEMENT  IV  LAVERS  ***** 

LAVER  STRESS  INTERVAL  INCREMENTAL 


NUOIER 

SETTLEMENT 

2 

1 

TO 

2 

-0.10946 

2 

2 

TO 

3 

0.21768 

2 

3 

TO 

4 

0.74912 

2 

4 

TO 

5 

9.5602S 

2 

5 

TO 

6 

-0.09385 

2 

LAVER  HISTORY 

1 .32375 

3 

1 

TO 

2 

-0.07111 

3 

2 

TO 

3 

0.09084 

3 

3 

TO 

4 

0.11109 

3 

4 

TO 

5 

0. 10559 

3 

5 

TO 

6 

-0.01898 

3 

LAYER  HISTORY 

0.21743 

5 

1 

TO 

2 

-0.032S4 

5 

2 

TO 

3 

0.04219 

5 

3 

TO 

4 

0.06297 

5 

4 

TO 

5 

0.07563 

5 

5 

TO 

6 

-0.01487 

5 

LAVER  HISTORY 

0.13338 

ttttt  settlement  iv 

STRESS 

INTERVAL  888ft 

STRESS  INTERVAL 

SETTLEMENT 

OVER  PROFILE 

1 

TO 

2 

-0.21311 

2 

TO 

3 

0.35070 

3 

TO 

4 

0.92318 

4 

TO 

5 

0.74147 

S 

TO 

6 

-0.12769 

TOTAL 

.  CLAV  SETTLEMENT 

1.67456 

(Continued) 
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Table  9  (Continued) 


t  CLAY  COMPRESSIBILITIES  < 


LAVER 

STRESS 

MU 

DELTA  E 

El 

E2 

REEDOUT 

15*  26  *  1$ 

2/13'79 

2 

1 

TO 

2 

0.00829 

-0.013S1 

0.85111 

0.86462 

2 

2 

TO 

3 

• .01259 

0.02686 

0.86462 

0.83776 

2 

3 

TO 

4 

0.02B5B 

0.09245 

0.83776 

0.74531 

2 

A 

TO 

S 

0.01801 

0.06914 

0.74531 

0.67617 

2 

S 

TO 

€ 

0.01570 

-0.01158 

0.67617 

0.68775 

3 

1 

TO 

2 

0.02145 

-0.01948 

0.78074 

0.80022 

3 

2 

TO 

3 

0.02085 

0.02489 

0.80022 

0.77533 

3 

3 

TO 

4 

0.01681 

0.03043 

0.77533 

0.74490 

3 

4 

TO 

S 

0.01300 

0.02893 

0.74490 

0.71597 

3 

S 

TO 

6 

0.01188 

-0.00520 

0.71S97 

0.72H7 

S 

1 

TO 

2 

0.01541 

-0.00267 

0.68394 

0.68661 

s 

2 

TO 

3 

0.01534 

0.00347 

0.68661 

0.68315 

s 

3 

TO 

4 

0.01492 

0.00517 

0.68315 

0.67797 

5 

4 

TO 

s 

0.01438 

0.00621 

0.67797 

0.67176 

5 

S 

TO 

6 

0.01418 

-0.00122 

0.67176 

0.67298 

t  SAND  SETTLEMENTS  I 

tsssstttttmstsssst 


tlttS  SETTLEMENT  IN  SAND  LAVERS  SSSSS 
THE  STRAIN  WEIGHTED  AVERAGE  BLOUCOUNT  -  33.3000 


LAVER 

NUMBER 

MEYERHOF 

4 

0.0090 

TOTAL  SAND 
SETTLEMENT 

0.0090 

TOTAL  CLAV 
SETTLEMENT 

1.6746 

TOTAL  PROFILE 

SETTLEMENT 

1.6835 

TIME-SETTLEMENT  RELATIONSHIPS 


IMICOEFFICIENT  OF  CONSOLIDATION*! If 
(CO) 


LAVER  NO. 
2 
3 
S 


CO 

SO. FT. 'DA V 
1.27M 

•.?••• 

0.4700 


SETTLEMENT  PER  LAVER  AT  THE  INPUT  INPUT 


TIMES 

SETTLEMENT 

DEGREE  OF  CONSILADTION 

DAYS) 

(FEET) 

LAVER  2 

(Uk) 

SO. 00 

-0.0331 

-2.  SO 

75.00 

0.0000 

0.00 

200.00 

0.4233 

31.98 

300.00 

0.1808 

66.  S4 

350.00 

1.0121 

76.45 

400.00 

1.1028 

83.31 

500.00 

1.2138 

81.69 

600.00 

1.2691 

95.87 

1000.00 

1.3804 

99.75 

(Continued) 


Table  9  (Concluded) 


REEDOUT 


16S26U5  2/13/79 


LAVER  3 

59.99 

<9.9366 

-16.84 

75.99 

<9.9173 

-7.95 

2««.N 

9.1931 

47.41 

3N.M 

9.2922 

93.99 

3SIM 

9.2144 

98. $9 

4N.M 

9.2163 

99.49 

5H.M 

9.2173 

99.93 

6H.M 

9.2174 

99.99 

1NMI 

9.2174 

199.99 

LAVER  S 

59.  •• 

<9.9988 

-6.57 

75.  M 

<9.9943 

-3.21 

20.99 

9.9394 

22.79 

399.99 

9.9772 

57.87 

359.99 

9.9997 

68.94 

499.99 

9.190 

75.58 

599.99 

9.1146 

85.95 

699.99 

9.1226 

91.92 

190.99 

9.1322 

99.12 

TOTAL 

SETTLEMENT  OF  PROFILE 

AT  TIMES  INPUT 

TIRES 

SETTLEMENT 

DEGREE  OF  CONSOLIDATION 

[  DAVS  ) 

(FEET) 

<U«) 

59.99 

-9.9784 

-4.68 

75.99 

-9.9215 

-1.29 

299.99 

9.5568 

33.25 

399.99 

1.1692 

69.28 

359.99 

1.3172 

78.66 

40.0 

1.4199 

84.79 

50.0 

1.5467 

92.31 

09.0 

1.6991 

96.0 

190.0 

1.670 

0.73 
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(Figure  7).*  Table  10  shows  the  results.  As  can  be  seen,  the  computer 
and  hand  solution  agree  very  closely. 

30.  Settlement  in  one  of  the  clay  layers  (layer  2)  determined 
from  MACSET1I  was  compared  to  hand  solutions  for  that  layer  using 
Terzaghi's  theory.  Table  11  shows  these  results.  The  computer  results 
compare  well  with  the  hand  solutions.  It  must  be  remembered  that  the 
theory  is  simplified  to  apply  to  the  soil  and  load  conditions,  so  these 
results  are  only  an  estimate. 

31.  Hand  calculation  of  the  degree  of  consolidation  of  t he  clay 
layers  was  based  on  the  methods  in  EM  1110-2-1904  for  clay  material. 

The  results  are  presented  in  Table  12.  There  are  no  noticeable  differ¬ 
ences  in  the  results  of  the  computer  and  hand  calculations. 

32.  Settlement  due  to  the  compression  of  the  sand  was  calculated 
by  Meyerhof's  method  using  a  vertical  strain  influence  factor.  The 
settlement  produced  by  MAGSETII  was  0.009  ft,  whereas  hand  calculations 
produced  0.007  ft,  values  which  are  very  close.  Again,  il  must  be  remem¬ 
bered  that  these  are  only  estimates  and  hand  calculations  require  some 
interpolation. 


Example  Problem  2 

33.  Example  Problem  2  was  taken  from  EM  1110-2-1904.  A  plan  view 
of  the  problem  is  shown  in  Figure  8.  Appendix  A  to  EM  1110-2-1904, 
which  describes  the  problem  and  includes  hand  computations,  is  included 
as  Figure  9. 

Input/output 

34.  Input  to  program  10016  is  shown  in  Table  13.  Output  from  the 
program  is  shown  in  Table  14.  Using  the  stress  information  provided  by 
program  10016,  input  to  program  MAGSETII  was  prepared  and  is  shown  in 
Table  15.  Output  from  this  program  is  included  in  Table  16. 


*  This  chart  is  from  notes  by  Prof.  Robert  D'Andrea  to  Course  No. 

CE3040,  "Soil  Mechanics,"  Worcester  Polytechnic  Institute,  Worcester, 
Mass. ,  1975. 
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SCALE  OF  DISTANCE,  OQ  =  DEPTH  Z 
AT  WHICH  STRESS  IS  COMPUTED 


gure  7.  Newmark's  influence  chart  for  vertical  pressure 
(influence  value  =  0.001) 
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Depth 

ft 


Table  10 

Comparison  of  10016  and  Hand  Solutions  for  Vertical 


Stresses  in  Example  Problem  1 


_ Vertical  Stress,  kips/ft _ 

10016  Solution  Hand  Solution 


i  0.882 

0.880 

0.509 

0.511 

0.103 

0.115 

Table  11 

Comparison  of  MAGSETII  and  Hand  Solutions  for  Settlement  in 


Layer  2  of  Example  Problem  1 


Increment 

No. 


_ _ Settlement, 

MAGSETII  Solution 

-0.10946 

0.21768 

0.74912 

0.56025 

-0.09385 


Hand  Solution 


Table  12 

Comparison  of  MAGSETII  and  Hand  Solutions  for  Rate  of  Settlement 


in  Layer  2  of  Example  Problem  1 


MAGSETII 

Solution 

Hand 

Solution 

Time 

Settlement 

Degree  of 
Consolidation 

Settlement 

Degree  of 
Consolidation 

days 

ft 

percent 

ft 

percent 

50 

-0.0331 

-2.5 

-0.037 

-3 

75 

0.000 

0.0 

0.000 

0 

200 

0.4233 

31.98 

0.50 

37 

300 

0.8808 

66.54 

0.98 

72 

350 

1.0121 

76.45 

1.08 

81 

400 

1.1028 

83.31 

1.21 

90 

500 

1.2138 

91.69 

1.25 

94 

600 

1.2691 

95.87 

1.29 

97 

1000 

1.32044 

99.75 

1.32 

99 

PART  (  XIX.  CHAPTER  4 
January  1953 


APPENDIX  A 

ILLUSTRATIVE  PROBLEM -SETTLEMENT  ANALYSIS 

The  Problem  Determine  the  total  settlement  resulting  from  a  buried  day  stiatum  and  the 
t imp-settlement  rate  for  a  structure  supported  by  nine  4-  bv  4-foot  footings  located  at  an  elocution 
It  feet  below  the  natural  ground  surface.  The  foundation  plan  for  the  structure  ,atid  the  sod  con¬ 
ditions  beneath  the  structure  are  shown  on  Plate  No  1.  The  gross  unit  load  on  each  fooling  is 
2  tons  per  sipiare  foot  The  construction  rate  of  load  w  ill  be  applied  uniformly  in  the  first  (it)  du\s 
and  the  remaining  25  percent  of  the  load  will  be  applied  uniformly  during  the  ne\t  :i0  dues  The 
amount  of  rebound,  resulting  from  the  excavation  before  the  construction  load  is  applied,  is  assumed 
to  be  negligible.  Consolidation  test  data  for  a  representative  sample  of  the  clut  stratum  are  shown 
on  Plates  Nos  2  and  K  on  width  are  plotted  pressure-void  ratio  ami  time-settlement  data,  respec- 
tivelv,  for  the  sample  tested  Examination  of  the  consolidation  data  shows  the  day  to  he  normallv 
consolidated 

Total  Settlement  In  order  to  determine  the  differential  settlement  to  be  expected  between 
footings,  the  settlement  must  be  computed  for  several  points  such  as  A,  B,  and  C,  Plate  No  I. 
Proceed  with  the  analysis  as  follows: 

(It  Construct  the  load-depth  diagram  for  the  existing  overburden  conditions,  Plate  No  i 
Since  the  clay  stratum  is  only  '20  feet  thick  it  is  safe  to  assume  that  the  pressure  distribution  i> 
uniform  from  top  to  bottom  and  that  the  pressure  at  the  middle  of  the  stratum  (depth  25  feet i 
represents  the  average  pressure  in  the  stratum.  Determine  ;n  from  the  load-depth  diagram.  Plate 
No.  3,  which  at  25  feet  is  1  15  tons  per  square  foot.  The  overburden  pressure,  pu  is  the  same  for 
all  three  points  A,  B,  and  C. 

(21  The  pressure  due  to  the.  added  load  of  the  structure  may  he  determined  either  as  point 
loads,  since  the  dimension-depth  ratio  is  greater  than  3.  or  as  area  loads.  The  area  load  method 
was  selected  for  this  problem  because  it  is  valid  for  ail  depths.  It  was  assumed  that  the.  Boussinesq 
solution  would  best  fit  the  soil  conditions  encountered  in  this  problem.  Plate.  No  0,  was  con¬ 
structed  ns  described  in  paragraph  4  03e  and  is  utilized  for  determining  the  vertical  pressures  at 
depths  of  10  and  25  feet  below  the  footings  for  each  of  the  three  points  A,  B,  and  C.  The  gross 
unit  load  of  the  footings  must  be  eorreeted  for  the  weight  of  5  feet  of  sand  which  was  excavated, 
to  obtain  the  net  load  applied  to  the  foundation.  From  Plate  No.  3  the  pressure  due  to  5  feet  of 
sand  is  0  31  ton  per  square  foot. 

Henee: 

5  =  2  00  —  0.31  =  1 .69  tons  per  square  foot  (net  load). 

Make  up  overlays  of  the  foundation  plan  shown  on  Plate  No.  1  with  scale  OQ  pqim)  to  10  and  25 
feet  for  use  with  Plate  No.  6.  Using  foundation  plan  with  chart  scale  equal  to  25  feet,  place  point 
A  over  the  center  of  the  chart.  Count  the  number  of  influence  areas  on  the  chart  which  fall  within 
the  outlines  of  the  individual  footing  areas.  Since  nil  footings  have  identical  unit  loads,  the  in¬ 
fluence  areas  under  all  footings  may  be  combined.  A  total  of  25  influences  is  counted.  Each 
influence  area  is  equal  to  0.001  q.  Then  for  point  A  at  a  depth  of  25  feet  below  the  footing  the 
pressure  due  to  the  structure  load  p,  is: 

p,=25X0.001X1.69=0.042  tons  per  square  foot. 


A-l 

Figure  9.  Appendix  A  to  EM  1110-2-1904  (Sheet  1  of  9) 
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Ym  ■  " 
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Yb  --  4  5#/Cu  ft 
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Figure  9.  (Sheet  2  of  9) 


PLATE  NO.  I 
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FAIT  cm.  CHAPTER  4 

Jaaaary  IMA 


Repeat  the  above  procedure  for  other  points  and  depths.  The  results  for  this  problem  are 
tabulated  as  follows: 


Depth  Mow  footing 

Point 

Number  of  influence 
areal 

P. 

tone  per  sq  ft. 

10 

A 

81 

0  136 

10 

B 

72 

.  122 

10 

C 

66 

112 

25 

A 

43 

.073 

25 

B 

33 

.  056 

25 

C 

25 

.  042 

Construct  a  load-depth  diagram,  Plate  No  6,  using  the  shove-computed  values  of  pressures. 

(3)  Using  Plate  No.  7,  determine  p,  at  the  middle  of  the  clay  stratum  which  is  located  20  feet 
below  the  bottom  of  the  footing  for  points  A,  R,  and  C.  Then,  p2  =  p,  1  p,.  After  obtaining  p2, 
determine  values  of  the  void  ratios  e,  and  r2  corresponding  to  p,  and  p2,  by  use  of  the  pressure-void 
ratio  curve,  Plate  No  2.  The  total  settlement  A II  may  be  obtained  by  the  relationship: 


AH= 


Cl  ~C2 
1  +«i 


2H, 


Results  of  the  foregoing  operations  are  summarized  below. 


Point 

A 

B 

c 

pi 

1.  15 

1.  15 

1.  15 

p 

.08 

.07 

.06 

pi 

1.23 

1.  22 

1.21 

<i 

1.044 

1. 044 

1.044 

1.035 

1.037 

1.038 

fi-ft 

.009 

.  007 

.  006 

l+e, 

2.  044 

2.  044 

2.  044 

2 11 

20 

20 

20 

AH 

.09 

.07 

! 

.06 

Differential  settlement  between  points  A  and  B  is  0.02  foot  and  between  points  A  and  C  is  0.03  foot. 

Time-tttUlemerU  ralt.  Using  Plate  No.  10,  the  time  for  50  percent  consolidation  of  the  labor¬ 
atory  specimen  1.25  inches  in  thickness  is  found  to  be  8.1  minutes  The  time  for  50  percent  con¬ 
solidation  of  the  field  stratum  is  then  found  by  the  use  of  the  relationship: 


Substituting: 

/  20  ft  \* 

(81  minutes) 

t„,=  298,598  minutes  =  207 .4  days. 


A-2 
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I'sing  the  relationship  in  equation  (IV)  as  follows: 


PAST  CX1X.  CHAPTER  * 
luurt  IMS 


t  '*7  T. 

r*  o.ioc 

compute  the  time  for  vurtous  percents  consolidation.  substituting  values  of  T„  found  in  Table  1. 
Substituting,  for  10  percent  consolidation 


MO/ 


298,598  min 
0.196 


>,0m077=  11  730  nun.  -n  14  days 


Results  of  the  computations  for  the  various  percents  d  consolidation  are  tabulated  below. 


f'% 

T. 

Minutes 

Dayt 

10 

0.  0077 

11,730 

8.  14 

20 

.  1)314 

47,  836 

32.  8 

30 

.  0707 

107,  70S 

74.  f 

40 

.  126 

191, 955 

123.  3 

50 

.  196 

298,  697 

207  3 

60 

.  286 

435,  709 

302  6 

70 

.  403 

613, 9o4 

426.  3 

80 

.  567 

863,801 

599.  9 

90 

.848 

I,  291,  893 

897.  1 

95 

1.  129 

1,  719,  985 

1194.  4 

Since  the  time  rate  for  construction  is  not  uniform,  it  should  be  divided  into  convenient  incre¬ 
ments  such  as  25  percent  of  the  load.  It  is  assumed  that  each  25  percent  increment  will  then  cause 
approximately  25  percent  of  the  total  consolidation,  and  the  settlement  resulting  from  each  incre¬ 
ment  is  assumed  to  start  at  the  half-time  for  the  application  of  increment  of  load  To  construct  a 
loading  diagram  as  shown  on  Plate  12,  divide  the  load  application  into  four  25  percent  increments. 
Draw  a  time-settlement  curve  for  each  increment,  starting  at  the  half-time  for  each  increment  and 
using  time  values  tabulated  above.  The  percents  consolidation  will  be  divided  by  4  in  each  case. 
The  time-settlement  curve  for  construction  loading  is  then  obtained  by  adding  the  ordinates  of  each 
increment  time-settlement  curve  and  drawing  a  smooth  curve  through  the  points  so  obtained. 
Adjust  the  initial  portion  of  the  curve  by  starting  the  curve  at  zero  time  and  settlement. 
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CORPS  OF  ENGINEERS 


Table  13 

Input  Data  File  for  Program  10016 


(Example  Problem  2) 


LIST  RLMBS 

1000  JAN  22, 1979  ) 

1010  SAMPLE  PROBLEM 

1020  9  4FT.  BV  4FT .  FOOTINGS  (  TtTLE 

1030  RLM 

1040  COMPARE  TO  EMU  10-2-1904  ' 

1060  1  9  (KODE.  NAREA) 

1070  1.00  0.0  0.0  0.0  0.0  4.0  4.0  4.0  4.0  0.0—-^ 
1080  1.00  0.0  18.0  0.0  18.0  4.0  22.0  4.0  22.0  0.0  \ 

1090  1.00  0.0  36.0  0.0  36.0  4.0  40.6  4.0  40.0  0.0 

1100  1.00  0.0  0.0  18.0  0.0  22.0  4.0  22.0  4.0  18.0 

1110  1.00  0.0  18.0  18.0  18.0  22.0  22.0  22.0  22.0  18.0 

1120  1.00  0.0  36.0  18.0  36.0  22.0  40.0  22.0  46.0  18.0 

1130  1.60  0.0  6.0  36.0  0.0  40.0  4.0  40.0  4.0  36.0 
1140  1.00  0.0  18  0  36.0  18.0  40.0  22.0  40.0  22.0  36.0 

1150  1.00  0.0  36.0  36.0  36.0  40.0  40.0  40.0  40.0  36.0 

1160  2  1  0.0  NDIST.  NWEST.  AMU 

1170  0.0  40.0  10.0  0.0  0.0  10.0  0.0  0 . 0 

1171  2  1  0  0 

1180  0.0  40.0  10.0  0.0  0.0  10.0  0.0  20.0  \ 

1181  2  1  0.0  \ 

1190  0.0  40.0  10.0  0.0  0.0  25.0  0.0  0  0  \ 

1191  2  1  0.0 

1200  0.0  40.0  10.0  0.0  0.0  25.0  0.0  20  0  / 

1210  0  0  0 . 0  —  (STOP  THE  PROBLEM)  / 


AINTL.  FINAL,  DELTA.  XP,  YP.  ZP,  SLP,  BLINE  - 

0(1),  ZLAY(1),  XC(1,1),  YC(1,1),  XC(2,1),  YC(2,1),  [ 
XC(3,1),  YC(3,1),  XC(4,1),  YC(4,1)  j 
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Table  14 

Output  Data  File  for  Program  10016 


RLMSB  1®-  3  59  2/14/79 

JAN  22. 1979 
SAMPLE  PROBLEM 
9  4FT  BY  4FT  FOOTINGS 
RIM 

COMPARE  TO  EMI  1 19-2-1984 
BOUSSINESQ  SOLUTION 

HORIZONTAL  STRESS  DISTRIBUTION  AT  DEPTH(Z)  •  19  89 


ELASTIC  SOLUTION  NORMAL  LOADING 


Y-COORDINATE 

x-coordinate 

UERTICAL  STRESS 

UERTICAL  STRESS 

e 

9 

8  864 

8  964 

8 

19  98 

8  837 

8  937 

8 

28  88 

9  872 

9  872 

8 

38  98 

8  837 

8  837 

8 

49  98 

8  864 

8  864 

NUMBER  OF  AREAS  USED  IN  CALCULATION  •  9 

NOTE-ALL  Z  VALUES  ARE  REFERENCED  TO  THE  LOWEST  PART  OF  THE  INPUT. 
CONFIGURATION 


JAN  22,1979 
SAMPLE  PROBLEM 
9  4 FT  BY  4FT  FOOTINGS 
RLM 

10MPARE  TO  EMU18-2-1994 
BOUSSINESQ  SOLUTION 


HORIZONTAL  STRESS  DISTRIBUTION  AT  DEPTH(Z)  •  18  89 


V-COORDINATE 

X-COORDINATE 

ELASTIC  SOLUTION 

VERTICAL  STRESS 

NORMAL  LOADING 

UERTICAL  STRESS 

29  99 

9 

9  972 

9  872 

29  99 

19  99 

9  942 

9  942 

29  99 

29  99 

9  982 

9  982 

29  99 

39  99 

9  942 

9  942 

29  99 

49  99 

9  972 

9  972 

NUMBER 

OF  AREAS  USED  IN 

CALCULATION  ■  S 

NOTE-ALL  Z  UALUES 

ARE  REFERENCED 

TO  THE  LOUEST  PART  OF 

THE  INPUT, 

CONFIGURATION 


(Cont inued) 


Table  14  (Concluded) 


jan  22.1979 

SAMPLE  PROBLEM 

9  4FT  BY  4FT  FOOTINGS 

RLM 

COMPARE  TO  EN1119-2-1994 
BOUSSINESO  SOLUTION 


HORIZONTAL  STRESS  DISTRIBUTION  AT  DEPTH(Z)  •  25  99 


Y-COORDINATE 

X-COORDINATE 

ELASTIC  SOLUTION 
UERTICAL  STRESS 

NORMAL  LOADING 
UERTICAL  STRESS 

9 

9 

9  923 

9  923 

9 

19  99 

9  927 

0  027 

0 

29  99 

9  929 

9  929 

9 

39  99 

9  92? 

9  927 

9 

49  99 

9  923 

9  923 

NUMBER  OF  AREAS  USED  IN  CALCULATION  •  9 

NOTE-ALL  2  UALUES  ARE  REFERENCED  TO  THE  LOUEST  PART  OF  THE  INPUT. 


RLMSB 


19 '  3  59  2/14/79 


CONFIGURATION 

JAN  22.1979 
SAMPLE  PROtLEP 
9  4FT  BY  4FT  FOOTINGS 
RLM 

COMPARE  TO  EH1119-2-1904 
BOUSSINESO  SOLUTION 

HORIZONTAL  STRESS  DISTRIBUTION  AT  DEPTH(Z)  ■  25  99 


ELASTIC  SOLUTION  NORMAL  LOADING 


Y-COORDINATE 

X-COORDINATE 

UERTICAL  STRESS 

UERTICAL  STRESS 

29  99 

9 

9  929 

9  929 

29  99 

19  99 

9  935 

9  935 

29  99 

29  99 

9  938 

9  938 

29  99 

39  99 

9  935 

9  935 

29  99 

49  99 

9  929 

9  929 

NUMBER  OF  AREAS  USED  IN  CALCULATION  •  9 

NOTE-ALL  Z  UALUES  ARE  REFERENCED  TO  THE  LOUEST  PART  OF  THE  INPUT. 
CONFIGURATION 
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Table  15 


Input  Data  File  for  Program  MAGS1CTI I 
(Example  Problem  2) 


SECOND  RUN 


THIRD  RUN 


Rinn 

9 '55  57 

1000 

1  2 

a  (PROBLEM  NO..  NO  OF  LAYER,  RUN  CONTROL) 

idle 

COflPARJSION  TO  £fU110 

-2-1904 

PROBLEN  (TITLE! 

1029 

1  1 

3  2 

1  i 

3  2 

(INPUT  OPTIONS) 

1630 

FEET 

TONS  (UNITSI 

1640 

0  0313 

15  0000 

(GROUND  WATER  DATA) 

1050 

1  0  0000  1 

1060 

2  IS. 

0000  1  (LAYER  INTERFACE  DATA) 

1070 

3  35.00000' 

1030 

1090 

1  3 

2  1 

0  0625  1 
0.0525  ! 

(SOIL  PROPERTIES) 

1109 

0 

0.  (DATUM  CONVERSION! 

1110 

1120 

1  2675 
0.4225 

®  |  (EFFECTIVE  STRESS  INCREMENT) 

1130 

2  0. 

0380  (EFFECTIVE  STRESS  DISTRIBUTION! 

1140 

2  1 

1.1500 

0  2400 

1159 

2  2 

1 . 1230 

0  5000 

1  1500 

0  2400 

0 

1160 

2  3 

1  0500 

1  0000 

1  1200 

0  5000 

0 

1170 

2  4 

6  9503 

2  0000 

1  0600 

1  0000 

0 

1180 

2  5 

0  8200 

4.0000 

0  9500 

2  0000 

1 

1190 

1 

(CONSOLIDATION  OPTION  CONTROL! 

1200 

2  095  0  0  (CONSOLIDATION  DATA) 

1210 

1  2 

0 

'1220 

COCIPRRISION  TO  EHllie 

1-2-1904 

PROBLEM 

1230 

1  1 

3  2 

1  0 

0  2 

1240 

FEET 

TONS 

1250 

0  0313 

IS  0000 

1260 

i  0  e 

1000 

1270 

2  IS 

0000 

1280 

3  35  00000 

1290 

1  3 

0  0625 

1300 

2  1 

0  0525 

1310 

0 

0 

1320 

1  2675 

0 

1330 

0  4225 

1 

1340 

2  0 

0290 

1350 

2  1 

1  1500 

0.2499 

1360 

2  2 

1.1200 

9  5099 

1  1500 

0  2400 

0 

1370 

2  3 

1.0690 

1  0000 

1.1200 

0  5003 

0 

1380 

2  4 

9.9500 

2  0000 

1.0600 

1  0000 

0 

1390 

2  5 

9  8209 

4  0000 

9.9500 

2  0000 

1 

1400 

1 

1410 

2  095  0  0 

1420 

1  2 

1 

"1430 

COrtPARISION  TO  EN1110-2-1904 

PROBLEfl 

1440 

1  1 

3  2 

1  0 

0  2 

1450 

FEET 

TONS 

1460 

0.0313 

15  0099 

1470 

1  0  0000 

1480 

2  15 

0000 

1490 

3  35  00000 

1500 

1  3 

0  0625 

1510 

2  1 

0  0525 

1520 

0 

0 

1530 

1  2675 

0 

1540 

0  4225 

1 

1550 

2  0 

9230 

1560 

2  1 

1  1599 

0  2490 

1570 

2  2 

1  1299 

0  5900 

1  1500 

0  2400 

0 

1580 

2  3 

1  0699 

1  0000 

1  1200 

0  5000 

0 

1590 

2  4 

0  9699 

2  0009 

1  0690 

1  3000 

0 

1600 

2  5 

9  8299 

4  0000 

9  9S99 

2  0030 

1 

1610 

1 

1620 

2  095  0  0 

2/14/79 


DEFORMATION 

CURVE 
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Table  16 

Output  Data  File  for  Program  MAGSETII 
(Example  Problem  2) 


RI.BW1  9.57.  g  2/14/79 


*  MAGSET-I1 

> 

*  MAGNITUDE  OF  SETTLEMENT  OF 
»  A  MULTl-LAVERED  SOIL  SYSTEM 

* 

t*tS**lttt*ttl*|It*||*t*t*I*tt 

iiiiiutiiimutinii 

*  SPECIFICATIONS  FOR* 

*  PROBLEM  NO  1  t 
Oil**  **********  tsisil 

*****  TITLE  ***** 
CONPARISION  TO  EMU1S-2-1BR4  PROBLEM 
*****  UNITS  ***** 

length  force 


FEET 

TON  S 

ttmmitiitm 

*  SOIL  PROFILE  t 

1  DESCRIPTION  « 

iiimuiiimxi 

DATUM  ELtUATION 

• 

DIFFERENCE 

IN  ELEVATION  • 

9 

LAVER  SOIL 

INTERFACE 

DATUM 

UNIT 

NUMBER  TYPE 

DEPTH 

elevations 

THICKNESS 

WEIGHT 

1  incqnp 

IS  •* 

t  tses 

IS  M 

-IB  M 

a  clav 

a*  •• 

•  9525 

35  •» 

-3S  *» 

UNIT  UEIQHT 

GROUND  WATER 

GROUND  WATER 

OF  WATER 

LEVEL 

DATUM  ELEVATION 

*  *313 

IS  •• 

'is  »• 

t  INSITU  EFFECT1UC 

tllttUIttlftttlll 

STRESS  t 

LAVER 

INPUT 

CALCULATED 

INSITU 

NUMBER 

VALUE 

VALUE 

STRESS 

1 

* 

»  4SIB 

e  hib 

2 

* 

l  149  5 

1  1495 

*»**i*tiiiti**tittt*t<t 

*  CLPV  SETTLEMENT  data  t 
*t*(*II*ltltt*lt *1*11*1* 

EFFECTIVE  STRESS  INCREMENTS  INPUT  IV 

A  DISTRIBUTION  FUNCTION 
*****  STRESS  INCREMENTS  FOR  STRATUM  ***** 

POINT  NUMBER  STRESS  INCREMENT 

1  i  as?s 

a  •  *aas 


*****  STRESS  DISTRIBUTION  FUNCTION  ***** 
LAVER 

NUMBER  VALUE 

a  *  tai* 


*****  EFFECTIVE  STRESS  HISTORV  1**11 
LAVE*  NO  PT  NO  STRESS  INCREMENT  STRESS  VALUES 


a  i 

a  a 


•  *412 

*  (1S1 


1  1*77 

I  2137 


(Continued) 
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Table  16  (Continued) 


RLNMM 


9 ■ 57 •  t  2/14/79 


DEFORMATION  CURVE6  INPUT  IV 

+  COORDINATE  POINTS 


mu 

COORDINATES 

OF  POINTS  ON  THE 

DEFORMATION  CURVES 

***** 

LAVER 

POINT 

REIOUND 

REBOUND 

NUMBER 

NUMIER 

VOID  RATIO 

STRESS 

VOID 

RATIO 

STRESS 

2 

1 

1  ISM 

9  2490 

2 

2 

1  12M 

0  S#00 

t 

1S00 

0  2400 

2 

3 

1  9600 

1  0009 

1 

1200 

0  5000 

2 

4 

0  9509 

2  0999 

1 

0600 

1  0990 

2 

5 

•  9200 

4  0000 

0 

9S90 

2  0000 

*****  SLOPES  ON 

the  deformation 

CURVES  ***** 

LAVER 

LINE 

cc 

CE 

CC 

CE 

NUMBER 

NUMBER 

(STRAIN) 

(STRAIN) 

2 

1 

0  9942 

0  0941 

0  0438 

0  0438 

2 

2 

0  1993 

0  2  993 

0  0940 

0  0940 

2 

3 

0  3654 

9  3654 

0  1774 

0  1774 

2 

4 

9  4319 

9  4319 

9  2215 

0  2215 

uiuiimutnuKitiummiii 

»  CLAV  SETTLEMENT  CONTRIIUTIONS  t 

nmiinmiiminnimmmi 


LAVER 

NUMIER 

2 

2 


*****  SETTLEMENT  »v  LAVERS  ***** 

STRESS  INTERVAL  INCREMENTAL 

SETTLEMENT 

1  TO  2  I  16393 

2  TO  3  I  12174 


I  LAVER  MISTORV  «  61467 


*****  SETTLEMENT  IV  STRESS  INTERVAL  ***** 


STRESS  INTERVAL 


TO 

TO 


TOTAL  CLAV  SETTLEMENT 


SETTLEMENT 
OVER  PROFILE 
•  16393 

I  12674 

I  11467 


***> **************** ****** 
«  CLAV  CONPRESSIIILITICS  * 
************************** 


LAVER 

STRESS 

MV 

DELTA  E 

El 

E2 

2 

1  TO  2 

0  06636 

•  S6ES1 

1  13769 

1  03138 

2 

2  TO  3 

9  99480 

6  6)211 

1  63131 

1  02926 

there  are  no  sand  lavers  in  the  soil  profile 


tttt(tttlttt|||tttl*tt!|ttlll!t|tt* 

•DEGREE  OF  CONSOLIDATION* 


••••COEFFICIENT  OF  CONSOLIDATION**** 

(CVS 


LAVER  NO 

I 


Cv 

SO  ft  /DAV 

I  6961 


(Continued) 


able  16  (Continued) 


9  57.  9 


DEGREE  OF  CONSOLIDATION 
U* 


TIflE  FACTOR 
7l> 


LAVER 


timtittiimimttit 

f  SPECIFICATIONS  FOR* 

*  PROBLEM  NO  1  * 

mtititmiitmttttk 


mttimtittm 
I  SOU  PROFILE  * 
*  DESCRIPTION  ( 

itiintmittm 


mu  TITLE  ***** 
COMPARISON  TO  EMllif-8-1904  PROBLEM 
*****  UNITS  ***** 

length  force 

FEET  TON  S 


8/14/79 

DAVS 


s 

99* 

•  0989 

8 

07 

19 

99* 

0  0070 

8 

88 

IS 

99* 

0  0177 

18 

63 

89 

99* 

0  0315 

33 

11 

85 

99* 

0  1498 

51 

74 

39 

99* 

0  0706 

74 

51 

35 

99* 

9  0963 

191 

41 

49 

99* 

0  1858 

138 

46 

45 

99* 

0  1593 

16? 

64 

59 

99* 

0  1966 

806 

07 

55 

00* 

0  8386 

851 

11 

69 

99* 

9  8863 

301 

35 

65 

00* 

0  3404 

358 

39 

79 

00* 

0  4028 

484 

05 

75 

00* 

6  4767 

501 

81 

89 

90* 

9  5671 

596 

99 

8S 

09* 

0  6837 

?19 

69 

99 

90* 

9  8480 

898 

63 

95 

99* 

1  1889 

1188 

87 

******* ttxttss*stttxttxs**ttxts 

*  * 

*  NAGSET-II  * 

*  X 

*  MAGNITUDE  OF  SETTLEMENT  OF  * 

*  A  flULTl-LAVERED  SOU  SYSTEM  * 

*  * 

«sttxtsxtx**t«ixxx**t ********** 


DATUM  ELEUATION 

0 

difference 

IN  ELEUATION  • 

0 

LAYER 

SOIL 

INTERFACE 

DATUM 

UNIT 

NUMBER 

TYPE 

DEPTH 

ELEUATIONS 

thickness 

WEIGHT 

0 

0 

1 

INCOflP 

15  90 

-15  00 

15  80 

0  0685 

8 

CLAY 

35  it 

-35  00 

80  00 

0  0585 

UNIT  UEIGHT 

GROUND  WATER 

GROUND  WATER 

OF  WATER 

IEUEL 

DATUM  ELEUATION 

9  0313 

IS  00 

-IS  00 

uuiiiiiiitiimitmuin 

*  INSITU 

EFFECTIUE 

STRESS  « 

itmtotiutusiiuutuii 

LAYER 

INPUT 

CALCULATED 

INSITU 

NUMBER 

UALUE 

UALUE 

STRESS 

1 

- 

0  4688 

0  4600 

8 

- 

1  1495 

1  1405 

(Cont inued) 
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Table  16  (Continued) 


RLMMM  9  57  8  2/14/79 

ts ttttt ttttttt ttttt ttttt 

*  CLAV  settlement  data  * 

t s sits ttttttt ttttt t ttttt 

effective  stress  increments  input  IV 

♦  A  DISTRIBUTION  FUNCTION 

Stitt  STRESS  INCREMENTS  FOR  STRATUM  ttttt 
POINT  NUMIER  STRESS  INCREMENT 

1  1  2675 

2  9  4225 


ttttt  STRESS  DISTRIBUTION  FUNCTION  ttttt 
LAVER 

NUMBER  UALUE 

2  9  9299 


ttttt  EFFECTIUE  STRESS  HISTORV  ttttt 
LAVER  NO  PT  NO  STRESS  INCREMENT  STRESS  UALUES 


2  1  9  9368  1  1863 

22  9  9123  1  1985 


DEFORMATION  CURUES  INPUT  fv 

COORDINATE  POINTS 


ttttt 

COORDINATES 

OF 

POINTS  ON  THE  DEFORMATION  CURVES 

ttttt 

LAVER 

POINT 

RESOUND 

REBOUND 

NUMBER 

NUMBER 

VOID  RATIO 

STRESS 

VOID  RATIO 

STRESS 

2 

1 

1  15*» 

t  2411 

2 

2 

1  12M 

•  Sitl 

1  1599 

9  240® 

2 

3 

1  tfM 

1  IMS 

1  1299 

9  5999 

2 

4 

•  «(«• 

2 

1  9699 

1  9999 

2 

5 

•  IMI 

4 

9  9599 

2  9999 

«»»•*  SCORES 

ON 

THE  DEFORNATION 

CURVES  ttttt 

LAVER 

LINE 

cc 

CE 

CC 

CE 

NUMBER 

NUMBER 

(STRAIN!  (STRAIN) 

2 

1 

•  *941 

t  SS41 

•  8438 

9  9438 

2 

2 

t  1993 

*  1993 

9  *94* 

B  9949 

2 

3 

•  3654 

•  3854 

•  1774 

9  1774 

2 

4 

•  4319 

(  4319 

S  221S 

9  2215 

t  ClAV  SETTLEMENT  CONTRIBUTIONS  t 

ttttttttttttttttttttttttttttttttt 


LAVER 

NUMBER 

2 

2 

2 


ttttt  SETTLEMENT  Bv  LAVERS  ttttt 

STRESS  INTERVAL  INCREMENTAL 

SETTLEMENT 

1  TO  2  9  94992 

2  TO  3  9  91699 

LAVER  HISTORV  9  96593 


ttttt  SETTLEMENT  BV  STRESS  INTERVAL  ttttt 
STRESS  INTERVAL  SETTLEMENT 

OVER  PROFILE 

1  TO  2  9  94992 

2  TO  3  9  91699 


TOTAL  CLAV  SETTLEMENT 


9  96693 


t  CLAV  COMPRESSIBILITIES  t 


(Continued) 
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Table  16  (Continued) 


rlnmh 


Layer 

STRESS 

WU 

2 

1  TO 

2 

1  166(1 

2 

2  TO 

3 

4  (654? 

DELTA  E 

*  99589 
«  88163 


SET  I  2/14/79 


El  E2 
1  *3799  1  *3289 
1  #3289  1  93128 


THERE  ARE  HO  SAND  LAVERS  IN  THE  SOIL  PROFILE 


iiitiiiunituitiiimmmuiiii 
TDESREE  OP  CONSOLIDATION* 
*******(******«************  ******* 


****COEFFICIENT  OF  CONSOLIDATION**** 
(CD) 


LAVER  NO 
2 


cv 

SO  FT  /DAV 

9  8958 


OESREE  OF  CONSOLIDATION  TINE  FACTOR  DAVS 

U*  TV 


LAVER  2 


S  99* 

8  8929 

2 

8? 

19  9m 

9  9979 

8 

28 

is  9m 

9  9177 

18 

63 

28  98* 

9  9315 

33 

11 

25  98* 

9  9492 

51 

74 

38  98* 

9  9788 

74 

Si 

35  99* 

9  9983 

181 

41 

48  99* 

9  1258 

132 

46 

45  99* 

9  1593 

1(7 

64 

5*  88* 

9  1988 

286 

97 

SS  98* 

9  2385 

251 

11 

89  99* 

9  2953 

381 

3S 

85  99* 

9  3484 

358 

38 

79  99* 

9  4921 

424 

85 

75  99* 

9  4787 

582 

81 

99  99* 

9  5871 

596 

99 

«  99* 

9  6837 

?19 

69 

84  99* 

9  8489 

892 

63 

95  99* 

1  1289 

1188 

27 

xxxxxxxxxxxxxxxxxxxixtxxtxxxxtx 

X 

X 

X 

NASSET-II 

X 

X 

X 

X 

HAGNITUDE  OF  SETTLENENT  OF 

X 

X 

A  HULTI-LAVERED  SOIL  SVSTEN 

X 

*  * 
********  t«  tt*X  t  **88**t  88*  **  *  I** 

*»****X»t*X*t(*tt*(* 

*  SPECIFICATIONS  FOR* 

*  PROILEN  NO  1  * 

*(*t**(***(t(tl*(t«**( 

*****  TITLE  ***** 
COHPARISION  TO  ENIH9-2-1994  PROILEN 


*****  UNITS  ***** 


X  SOIL  PROFILE  X 

X  DESCRIPTION  I 

IENCTH 

FEET  TON  S 

FORCE 

DATUN  ELEVATION  • 

8 

DIFFERENCE  IN  ELEVATION  • 

9 

(Continued) 
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Table  16  (Continued) 


RLHMH 

9  57  8 

2/14/79 

layer  soil 

INTERPACE 

DATUH 

UNIT 

NUMBER  TYPE 

depth 

ELEVATIONS 

thickness 

WEIGHT 

• 

• 

I  INCOHP 

IS  M 

-15  99 

15  99 

9  6625 

a  CLAY 

35  ee 

-35  19 

20  09 

9  0525 

UNIT  WEIGHT 

GROUND  WATER 

GROUND  WATER 

OP  WATER 

LEVEL 

DATUH  ELEVATION 

0  0313 

15  00 

-IS  90 

********** ******* ********** 

*  INSITu  EPFECTIUE 

STRESS  I 

miimiitiiimtnimm 

LAYER 

INPUT 

CALCULATED 

INSITU 

NUMBER 

UALUE 

UALUE 

STRESS 

l 

- 

0  4688 

0  4688 

2 

“ 

1  1495 

1  1495 

timnmmuiiimitt 

*  CLAY  SETTLEMENT 

DATA  * 

ittitmnmitittumt 

EFFECTIVE  STRESS  INCREMENTS  INPUT  IV 

A  DISTRIlUTION  FUNCTION 

i»«*«  stress  increments  for  stratum  t>>» 

POINT  NUMIER  STRESS  INCREMENT 

1  1  8675 

2  I  422S 

*»«««  STRESS  DISTRIlUTION  FUNCTION  !l!» 

LAVER 

NUMIER  VALUE 

2  •  R23R 


Hill  EFFECTIUE  STRESS  MISTORV  IKK 
LAVER  NO  PT  NO  STRESS  INCREMENT  STRESS  VALUES 

2  1  2  12*2  I  J7S7 

2  2  *  1**7  1  18*4 


DEFORMATION  CURVES  INPUT  |V 

COORDINATE  POINTS 


ttm 

COORDINATES 

OP 

POINTS  ON  THE  DEFORMATION  CURVES 

***** 

LAYER 

POINT 

REBOUND 

REBOUND 

NUMBER 

NUMBER 

VOID  RATIO 

STRESS 

UOID  RATIO 

STRESS 

2 

1 

1  1500 

0  2400 

a 

a 

1  1200 

0  5000 

1  1500 

0  2400 

a 

3 

1  0600 

1  0000 

1  1200 

0  5000 

a 

4 

0  9500 

a  0000 

1  0600 

1  0000 

a 

5 

0  8200 

4  0000 

0  9S00 

a  0000 

*****  SLOPES 

On 

the  depormation 

CURUES  IXttl 

layer 

LINE 

cc 

CE 

c.c 

r.F 

NUMBER 

NUMBER 

(STRAIN)  (STRAIN) 

a 

1 

0  0941 

0  0941 

0  0438 

0  0438 

a 

2 

0  1993 

0  1993 

0  0940 

0  0940 

a 

3 

0  3654 

0  3654 

6  1774 

0  1774 

2 

4 

0  4319 

0  4319 

9  2215 

0  2215 

(Continued) 
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Table  16  (Concluded) 


RLMMN 


9  57  8  2/14/79 


itmtiimiiiiimiiiniisiimi 

»  CLAV  SETTLEMENT  C0NTRX1UTI0NS  * 

tSXXlXSXXIXtltXXHXXXXtXXXXtXXXXX 


LAVER 

NUNSER 

2 

2 

2 


MX**  SETTLEMENT  IV  LAVERS  ***** 

STRESS  INTERVAL  INCREMENTAL 

SETTLEMENT 

1  TO  2  I  63911 

2  TO  3  0  61279 

LAVER  HISTORV  6  65179 


*****  SETTLEMENT  IV  STRESS  INTERVAL  ttttt 


STRESS  INTERVAL 


TO 

TO 


TOTAL  CLAV  SETTLEMENT 


SETTLEMENT 
OVER  PROFILE 
0  63961 
6  61279 

6  65179 


mmmtttmmtiiimi 
*  CLAY  COMPRESSIBILITIES  * 

iiimitittiiiimiumix 


LAVER 

STRESS 

(1U 

DELTA  E 

El 

E2 

2 

1  TO 

2 

9  96681 

9  99397 

1  03789 

1  03391 

2 

2  TO 

3 

9  96593 

9  99139 

1  63391 

1  03261 

THERE  ARE  NO  SAND  LAVERS  IN  THE 


SOIL  PROFILE 


mmtimittmmmunmtmt 

(DECREE  OF  CONSOLIDATION* 

XXXtX(XStt(XX*X*X*X***ttSSX*SX*tt(C 


XXXXCOEFFICIENT  OF  CONSOLIDATION**** 

(CV> 

CV 

LAYER  NO  SO  FT  /DAY 

2  6  6956 

DECREE  OF  CONSOLIDATION  TIME  FACTOR  DAYS 

UR  TV 


LAVER  2 


5 

00* 

9  9939 

2 

07 

10 

00* 

9  9979 

8 

28 

15 

00* 

9  9177 

18 

63 

20 

00* 

9  9315 

33 

U 

25 

00* 

9  9993 

51 

74 

30 

00* 

9  9798 

74 

51 

35 

06* 

9  9963 

101 

41 

40 

00* 

9  1358 

132 

46 

45 

00* 

9  1593 

16? 

64 

50 

00* 

9  1966 

206 

97 

55 

06* 

9  3386 

251 

11 

60 

66* 

9  3863 

361 

35 

65 

66* 

9  3999 

358 

30 

70 

00* 

9  9938 

424 

05 

75 

60* 

9  9767 

501 

81 

86 

00* 

9  5671 

596 

99 

86 

06* 

•  6837 

719 

69 

96 

00* 

8  8989 

892 

63 

95 

00* 

1  1889 

1188 

27 
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Comparison  of  results 

35.  The  stress  induction  and  settlement  results  for  example 
problem  2  from  EM  1110-2-1904  are  compared  with  the  computer  solutions 
in  Tables  17  and  18.  It  appears  that  although  EM  1110-2-1904  states 
that  the  values  are  calculated  at  25  ft  below  the  footing,  they  really 
seem  to  be  computed  at  20  ft  below  the  footing.  The  comparisons  are 
very  close. 

36.  The  time-settlement  results  are  compared  in  Table  19.  Again, 
the  results  are  very  close  showing  the  validity  of  the  programs  used. 

Table  17 

Comparison  of  10016  and  EM  1110-2-1904  Solutions  for 
Stress  Induction  in  Example  Problem  2 


Stress  Induction 

10016 

Solution  „ 

(Output  x 

1.69  tons/ft  ) 

EM 

1110-2-1904 

Location 

@  25  ft 

@  20  ft 

Solution 

A 

0.064 

0.073 

0.073 

B 

0.049 

0.058 

0.056 

C 

0.039 

0.046 

0.042 

Table  18 

Comparison  of  MAGSETII 

and  EM  1110-2-1904  Solutions 

for 

Settlement 

in  Example  Problem  2 

Settlement , 

ft 

MAGSETII 

Solution 

EM  1110-2-1904 

Location 

<?  25  ft 

@  20  ft 

Solut ion 

A 

0.085 

0.094 

0.09 

B 

0.065 

0.071 

0.07 

C 

0.052 

0.063 

0.06 

Table  19 

Comparison  of  MAGSETII  and  EM  1110-2-1904  Solutions  for  Rate 
of  Settlement  in  Example  Problem  2 


Degree  of 
Consolidation 


MAGSETII 

Solution 


Time,  days 


EM  1110-2-1904 
Solution 


18.63 

33.11 

51.74 

74.51 

101.41 

132.46 

167.64 

206.97 

251.11 

301.35 

358.30 

424.05 

501.81 

596.99 

719.69 

892.63 

1188.27 


133.3 


207.3 


302.6 


426.3 


599.9 


897.1 

1194.4 
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PART  IV:  EXAMPLE  PROBLEM  ILLUSTRATING 
INPUT/OUTPUT  FOR  PROGRAM  FD31 


37.  This  Part  contains  an  example  problem  solved  using  Program 
FD31.  Other  examples  solved  using  FD31  are  presented  in  Olson. 

Input 

3d.  Data  input  falls  into  three  categories:  (a)  input/output 
control  parameters,  (b)  raw  data,  and  (c)  program  control.  The  input/ 
output  control  parameters  govern  the  form  in  which  data  are  input  and 
the  amount  of  output  and  type  of  output  to  be  printed  out.  Raw  data 
consist  of  water  table,  drainage  data,  effective  weights,  layer  depth, 
excess  pore  pressures,  embankment  description  or  load  description,  time 
table  of  construction  for  detail  output,  consolidation  curve,  and  coef¬ 
ficient  of  permeability  and  coefficient  of  consolidation.  The  program 
control  data  consist  of  parameters  controlling  the  accuracy  of  computa¬ 
tions  within  a  program. 

39.  The  example  problem  soil  profile  (Figure  10)  is  a  simple  one 
for  which  Terzaghi's  theory  is  applicable.  The.  soil  system  consists  of 
1  ft  of  incompressible  sand  over  10  ft  of  clay  over  incompressible  sand 


APPLIEO  FILL 


•  -SAND  ' 

CLAY 

— ?  7T  - 7T - r  —  —  —  r  7V - 7'  -  7K — T - 

y'\/\  V"-..-''' 


Figure  10.  Soil  profile  for  FD31  example  problem 


All  layers  are  saturated  and  have  submerged  unit  weights  of  50  lb/ftJ. 

The  clay  is  linearly  elastic  and  the  consolidation  curve  passes  through 

the  points  (e,o)  =  (2.0075,  0  lb/ft2)  and  (1.9500,  2300  lb/ft2).  The 

2 

clay  has  a  constant  coefficient  of  consolidation  of  0.05  ft  /day).  Con¬ 
solidation  results  from  the  application  of  20  ft  of  fill  at  time  zero. 

2 

The  fill  has  a  total  unit  weight  of  112.4  lb/ft  and  is  also  saturated. 

The  water  table  is  30  ft  above  the  original  ground  surface. 

40.  The  complete  input  is  shown  in  conversational  mode  in  Table  20. 
Values  of  various  control  parameters  are  indicated  below  with  brief 


explanations : 

a .  101-109.  All  10  parameters  are  set  to  1  for  this  example 

in  order  to  obtain  maximum  output  information. 

Id .  K0DWT=C.  The  water  table  will  be  maintained  at  a  constant 
elevation. 

£.  T0PB=F.  The  upper  boundary  of  natural  soil  is  freely 
draining. 

d.  B0TB=F.  The  bottom  boundary  is  also  freely  draining. 

e.  USEQ0=T.  The  value  of  this  parameter  is  not  relevant 
for  this  problem  and  may  be  set  at  either  T  or  F. 

_f .  SANDPP=F.  The  value  of  this  parameter  is  not  relevant 
either  but  should  be  set  equal  to  either  F  or  S . 

£.  ALMX1=Q. 5.  The  value  of  this  parameter  is  irrelevant  to 
this  problem  because  the  loading  is  a  single  step  loading. 

]i.  ALMX2=0. 5.  This  is  the  upper  limit  on  all  values  of  A(IN) 
at  time  zero. 

i.  ALMX3=20. 5.  This  is  the  upper  limit  on  values  of  A(IN)  at 
time  TL(JLFIN);  i.e.,  the  last  point  on  the  loading  curve. 

J^.  CHGMIN^O. 1 .  The  value  of  CHGMIN  is  irrelevant  because  the 
coefficient  of  consolidation  is  constant. 

k.  CHGLIM=1 . 0.  The  value  of  this  parameter  is  irrelevant 
because  of  the  constant  coefficient  of  consolidation. 

l.  FCV=1 .0.  The  value  of  this  parameter  is  also  irrelevant 
because  of  the  constant  coefficient  of  consolidation. 

m.  TOI,=0. 1 .  The  Gauss-Seidel  subroutine  will  continue  iterat¬ 
ing  until  no  excess  pore  pressure  in  the  system  changes  by 
more  than  0.1  lb/ft^  on  the  last  iteration. 

n.  ITERMX=100.  If  the  number  of  iterations  in  the  Gauss- 
Seidel  subroutine  reaches  100,  the  calculations  will  be 
aborted  and  the  analysis  stops. 
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Table  20 

Input  In  the  Conversational  Mode  for  Program  FD31  Example  Problem 


IF  DATA  IS  TO  EE  PERU  FROM  A  DATA  FILE 
It  (PUT  THE  FILE  HOME  <8  C HARRCTERS  MAX.) 

HIT  CARRIAGE  RETURN  IF  INPUT  IS  FPDM  TERMINAL 


IF  DATA  IS  TO  EE  SAVED  TO  A  DATA  FILE 
INPUT  THE  FILE  NAME  <8  CHARACTERS  MAX.) 

HIT  CARRIAGE  RETURN  IF  NO  FILE  IS  TO  EE  WRITTEN. 


IF  OUTPUT  IS  TO  EE  WRITTEN  TO  A  FILE 
INPUT  THE  FILE  NAME  <8  CHARACTERS  MAX.) 

HIT  CARRIAGE  RETURN  IF  OUTPUT  IS  TD  COME  TD  TERMINAL. 

DO  YOU  WISH  TO  HAVE  ALL  INPUT  ECHOED  ?  (Y,'N> 

=N 

INPUT  TITLE  FOP  THIS  PUN 
=EX AMPLE  PROBLEM  FOR  SETTLEMENT  PACT  AGE 
ASSIGN  VALUES  OF  1  IF  OUTPUT  IS  DESIRED.  OTHERWISE  0. 

101  =  RESIDUAL  POPE  PRESSURES  AND  EFFECTIVE  STRESS 

102  =  LOAD  HISTORY 
IDS:  =  OUTPUT  TIMES 

104  =  E-P  CURVES 

105  =  INITIAL  VALUES  OF  CV 

1 08  =  0  FOR  DATA  AT  MIDHEIGHT  OF  LAYERS  AND  AT  INTERFACES 
BETWEEN  LAYERS 

=  1  FOR  DATA  AT  ALL  NODES  AND  INTERNDDES 
107  =  1  IF  DATA  APE  TO  EE  OUTPUT  ONLY  AT  TIMES  TD  •' JD> 

=  0  IF  DATA  ARE  TD  EE  OUTPUT  AFTER  EVERY  TIME  TC .  IN  THE 
CASE.  106  IS  AUTOMATICALLY  SET  AT  1 
IDS  =  0  FOR  NO  OUTPUT  OF  CONTROL  PARAMETERS 
=  1  FOR  OUTPUT  OF  CONTROL  PARAMETERS 
101. 102. I  03 .104.  1 05 . I  06 .107,  1 08 
=11  1  11111 

ZWTO=DEPTH  FROM  ORIGINAL  GROUND  SURFACE  TO  THE  WATER  TABLE. 
POSITIVE  DOWNWARDS  <FEET) 

ZWTF=DEPTH  FPOM  A  DATUM  ESTABLISHED  AT  THE  ELEVATION  OF  THE 
ORIGINAL  GROUND  SURFACE  TO  THE  FINAL  WATER  TABLE, 

POSITI VE  DOWNWARDS  (FEET) 

KQDWT=G  MEANS  THE  ELEVATION  OF  THE  WATER  TABLE  IS  CONSTANT 

R  MEANS  THE  ELEVATION  OF  THE  WATER  TABLE  RISES  FROM 

EWTO  TO  EWTF  AS  A  RESULT  OF  CONSTRUCTION  OPERATIONS 

D  MEANS  THE  ELEVATIDN  OF  THE  WATER  TABLE  DROPS  FROM 

EWTO  TO  EWTF  INSTANTLY  AT  TIME  ZERO 
S  MEANS  THE  ELEVATION  OF  THE  WATER  TABLE  ALWAYS  IS  THE 
ELEVATION  OF  THE  ORIGINAL  GROUND  SURFACE  AS  THAT 
SURFACE  SETTLES  OP  PISES. 

TORE  =  I  FOP  AN  IMPERVIOUS  UPPER  BOUNDARY 

=  F  FOP  A  FREELY  DRAINING  UPPER  BOUNDARY 
BOTB  =  I  FOP  AN  IMPERVIOUS  LOWER  BOUNDARY 

=  F  FOR  A  FREELY  DRAINING  LOWER  BOUNDARY 
USEQ0=  T  IF  STATIC  PORE  PRESSURES  ARE  ZERO  ABOVE  THE  WATER 
TREBLE  AT  ALL  TIMES 

*  F  IF  STATF'P  *  I)  =  1  Z  <  I )  -ZWT  >  ♦62. 4  AT  ALL  TIMES 
SAN DPP  =  F  IF  ALL  SAND  LAYERS  DRAIN  FREELY  HOP  I ZONT ALLY 
=  S  IF  SAND  LAYERS  APE  SEALFD  HORIZONTALLY 
ZWTD , ZWTF , KODWT , TC  PB , TOPB , BOPB»  U  S EC  0 . S ANDPP 
=-30.0  -30.0  C  F  F  T  F 


(Continued) 


Table  20  (Continued) 


HL=NUNBER  OF  LAYERS  OF  SOIL  'HOT  C  DUHT I  HR  MEM  FILL:' 

HL'L'=THIO  MET:  OF  LAVER  L  RT  TIME  ZERO  'FEET> 

GAMF'P 1  L>  =SUB:MERGED  UNIT  WEIGHT  OF  "OIL  IN  LHVEF'  L  HT  TIME  ZEF'Q 
•  PC  F  • 

MZ  'L'  =  NUMBEP  OF  NODES  IN  AMY  LHVEF'. 

NUMBER  OF  LHYER 

INPUT  HL  <  GAMFR  *  MZ  FOP  LHYEP  £ 

=  1.0  0  5  0 . 0  £ 

INPUT  HL ' GANPR. MZ  FOP  LHYEP  £ 

=10. 0  5 0 .0  11 

ALMX1  =  MH: :  I  MUM  HLPHH  HUP  IMG  LOADING  OP  UNLOHIiING 
ALMX2  =  MAXIMUM  HLPHH  HUP  I MG  A  NON-CONSTRUCTION  PEP  I OH  EXCEPT 
FOP  TC.  GT.TL  •:  JLEIM-1> 

HLMX3  =  MAXIMUM  HLPHH  HT  TC-TL • JLF IN • 

CHGLIM  =  MAXIMUM  ALLOWABLE  CHANGE  IN  CV  EXP PETS EH  HI  H  PAT  ID 
CHGM1N  =  THE  PPDGPAM  M ILL  NOT  CALCULATE  AVERAGE  VALUE'  OF 
CV  AND  RECYCLE  IF  CHGMHX. LT. CHGMIM.  RELEVANT  ONLY  FOP  A 
VARIABLE  PROPERTIES  SOLUTION. 

TOL  =  MAXIMUM  ALLOWABLE  CHANGE  IN  POPE  PRESSURE  BETWEEN 
ITERATIONS. 

I  TERM X  =  MAXIMUM  NUMBER  OF  ITERATIONS  IN  SOLVE 
HTML  IN  =  UPPER  LIMIT  ON  TM  TIME  DURING  A  TIME  HT 
FC V=A  REDUCTION  FACTOR  USED  WHEN  CHGMHX. GT. CHGLIM 
C H G P M X = M A X I M U M  ALLOWABLE  CHANGE  IN  EFFECTIVE  STRESS  IN 
SUBROUTINE  CONSOL  WITHOUT  ITERATING  FURTHER 

PPL IM=A  NUMBER  SUCH  THAT  THE  PROGRAM  QUITS  RUNNING  WHEN  ALL 
VALUES  OF  PP  APE  LESS  THAM  THIS  VALUE 
ALMX1 .  ALMX2?  AMLX3 

—  C*  B*  2  |"|  |*| 

chgmin.’chglIm.fcv 

= .  1  1.0  1.0 

TOL  . I  TER MX . HTML  I M - CHGMX -  PPL  I M 

= .  1  100  .  0  .1  .1 

KODPPP  =  C  IF  RESIDUAL  PORE  PRESSURES  APE  ALL  THE  SAME 
=  V  IF  RESIDUAL  PORE  PRESSURES  APE  INPUT 
AT  EACH  DEPTH  NODE 

LOBF'PP 

PPEPPjR.ES  I  DUAL  POPE  PRESSURE 

=  0.  0 

NFL  =  NUMBER  OF  CONSTRUCTION  STAGES .  MAY  BE  ZERO 
TFC BEG  < JF >  =T I ME  FILL  CONSTRUCTION  BEGINS  .  DAYS 
T  F  C  END  J  F  >  =  T I  ME  FILL  CONSTRUCTION  ENDS.  DAYS 
TL ' JL •  =  T I ME  WHEN  FILL  ELEVATION  IS  DEFINED.  DAYS 
EFL' JL ■  =ELEVATION  OF  THE  TOP  OF  THE  FILL  RELATIVE  TO  THE  TDP 
OF  THE  ORIGINAL  GROUND  SURFACE 
GF  <  JF>  =  TOTAL  UNIT  WEIGHT  OF  FILL.  F'C  F  . 

NFL 
=  1 

TFC BEG. TFC END. GF  1 
=  0.  00  0.  1  12.4 

TL «  EFL  1 
=  0 .  0  0 .  0 
TL «  EFL  £ 

=  0.  0  £0.  0 
TL • EFL  3 
*5000. 0  £0. 0 
TL.EFL  4 

=-l . 0  0. 0  (Continued) 
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Table  20  (Continued) 


TO  ■  JO  ■  =TIMES  FOP  MHICH  OUTPUT  IS  DESIRED 
I  OUT  •:  JO  •  =B  FOP  BETRILEB  OUTPUT  RT  TO  •:  JO  •  -M  FOP  HOME 
OTDRTH=B  FOP  DU T PUT  OF  DETAILED  BRTR  RT  THE  LR8T  VRLUE  OF  TC 
WHEN  RLL  PP ■  I  '  RPE  LESS  THRU  PPL  I M 
TO. I OUT  1 
=  .  1  D 
TO. I  OUT  £ 

=  1.0  D 
TO. I OUT  3 
=4.  0  M 
TO.  I  OUT  4 
=  1  0 .  D 

TO. I OUT  5 
=40.0  M 
TO.  I  OUT  6 

=  1  0  o .  o  D 

TO. I  OUT  7 
=4  0  0.  0  N 
TO. I  OUT  8 
=  1  0  0  0 .  0  I1 
TO. I  OUT  3 
=8000.0  D 
TO.  I  OUT  1  0 
=-  1  .  0  ft 

OTBRTR 
=  D 

FOBOPB=V  FOP  VO IB  PRTIO  D I RGPRMS 
=S  FOP  STRAIN  B I RGPRMS 

FOB  I NT  =  H  IHTEPPOLRTE  P  RP I THMET I C RLL Y  FDP  E-P  PELRT  I  DM'"  HI  P 
=  L  IHTEPPOLRTE  F  LOGRPITHMICRLLY  FDP  E-P  PELRT I  OH  I  HIP 
I OB0PB« I OBIMT 
=V  H 
LHVEP  1 

PE  =EFFEC T I VE  STRESS. EF=VOIB  PRTIO  OP  STRAIN. PIOHT  1 
=  0 .  0  0  1 .  0  0  O 

PE=EFFECTIVE  STPESS. EP=VO IB  PRTIO  OP  STPRIH.PIOHT  8 

=  1  0  0  0  0 .  0  0  1 . 0  0  0  0 

pe Effective  stress. ep=vdib  prtio  dp  strain. pi ont 

=  -  1  .  o  0.0 

i:cp=:lope  of  peboijhb  curve 
= ,  0  0 1  0  0 
LR'i'EP  8 

PE=EFFEC T IVE  STRESS. EP=VOIB  PRTIO  OP  STPRIH.PIOHT  1 

=  0  .  0  0  8 . 0  075 

PE  =  EFFEC  T  IVE  STRESS .  EF'=VO  I B  PRTIO  DP  STPRIH.PIOHT  8 
=8300.0  1.9500 

PE=EFFEC T IVE  STRESS > EP=VO ID  PRTIO  OP  STPRIH.PIOHT  3 
=  -  1 . 0  i'i  0 .  0  0 

CCP=:lOPE  OF  PEBOUHB  CURVE 

= .  0  0 1  u  0 

cv:b=r  LIMIT  SUCH  THhT  RHV  LR-.EP  mi TH  R  HIGHER  VRU IE  DF  cv 
IHOULB  BE  R  SRHB  LRYEP. 

c  v :  b 
=  1  o  o .  o 
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Table  20  (Concluded) 


t  ODSF'  ■  L>  =  c  FOP'  CON!  TAM T  CV  HMD  PK  VALUE  I 

=  V  FOR  CV  AMD  Ft  FUNCTIONS  OF  EFFECTIVE  STPEM 
VP  =  CONST  IF  ALL  LAYERS  HAVE  t QDSP vL' =C 
=  VAR  IF  AMY  LAYER  HAS  t ODSR • L> =V 
PCVCL- JCV> ^EFFECTIVE  STRESS  HT  WHICH  A  POIMT  DM  THE  R-CV 

CURVE  IS  DEFINED  FOR  LAYER  L  AMD  FOIMT  JCV  'RIF- 
C  VP  •  L  «  JC  V  >  =C  DEFF I C I ENT  OF  CONSOLIDATION  AT  THE  JCV  POIMT 
ON  THE  INPUT  F'-CV  CURVE  OF  LAYER  L  •  SO.  FT.  DAY 1 
Ft  P >  L ’  JC  V>  =C OEF  F I C I ENT  OF  PERMEABILITY  AT  THE  JCV  POINT  DM  THE 
INPUT  P- PR  CURVE  OF  LAYER  L  iFT.  DAY' 

LAYER  1 

C  VP  .  Ft  P 
=£00.0  £00. 0 
LAYER  £ 

CVP.PFP 
=  .  05  . 05 


lng ut  Data  Fi le  for  Program  FIVJl 
(Example  Problem  1) 


1  Oft 

:mHFLE  FFQtLEM  r-G. 

1  K« 

1  1  1  l  1 

2  2  l 

Id  ft 

- }  0 .  o 

-  *.  ft .  ft 

1  ?  ft 

i- 

140 

1 .  0  ft 

«.  ft .  ft 

i  *;  o 

1  0  .  0  0 

e.  ft.  ft 

i  *••  o 

e;  cj 

£0.  0 

i  ro 

.10  1  .  »'»!> 

1 . 0 

1.5  0 

.1  10  fi 

.ft  .1 

l  =*“ 

c 

c '  1 0 

ft.  ft 

f  o 

1 

CC  0 

0.  00 

0.  ftft 

c  ?.■:» 

0.  0  ft 

ft.  ftft 

c  4  0 

0.  ft  0 

j  ft.  ftft 

c.  i’t 

*.  ft  ft  ft.  00 

eft.  ftft 

c-  0 

-  1 .  ft  ft 

-ft.  ftft 

.*-  ~  ft 

.  i  r> 

„  -  0 

i .  o  p 

L  •-  I) 

4.  0 

**u 

lft.0  I' 

:  i  o 

4  IT  .  ft 

:  fit 

I  ft  ft  .  ft  fl 

••  •  IT 

4  ft  ft  ,  ft 

:'4  0 

1  ft  i»  ft .  ft  p 

0 

c  ft  ft  ft .  0  P 

;•»*.  »T 

-1.0 

'  ~  II 

D 

V  Cl 

h 

'  ift 

ft,  no 

1  ,  ft  ft  ft  ft 

4  0  0 

1  ft  ft  ft  ft ,  ft  ft 

1 .  111? II ft 

4  l  0 

- 1 .  ft  o 

-2 .  ftOftf 

4c- 0 

•  ft  ft  1  Oft 

4  SO 

ft.  Om 

£.  oor?. 

44  0 

*  Soft,  i  ni 

! .  A^Oft 

4**0 

-1 .  ft  ft 

-ft.  ftft flfl 

4*.  0 

.  0  ft  1  ft  ft 

4T0 

1 .  ft  ftE  ♦  Oc 

45?  0 

c 

4-40 

c.  ft  ftE  ♦ft*?  c . 

ftllg  ♦  il£ 

*.on 

c 

*:  t  it 

ft  ftp- ft  c  *. 

ft  llF  -  Oe' 

fvu  o:  tr<i?  f  i* :  i  ■ 

c  f  .*  f  r  f 

i  i 

i 

I  1  £ .  4  0 


73 


£.  DTMLIM=5 . 0.  The  analysis  will  be  aborted  if  the  running 
time  for  any  one  value  of  time,  TC,  reaches  5  seconds. 

£.  CHGPMX=0. 1 .  The  equation  for  effective  stress  at  any 
node  requires  as  input  the  value  of  surface  settlement, 
but  this  is  the  settlement  to  be  calculated.  The  sub¬ 
routine  will  iterate  on  the  effective  stresses  until  no 
effective  stress  changes  by  more  than  0.1  lb/ft^  during 
the  final  iteration. 

£.  PPLIM=0. 1 .  During  the  time  period  between  TL(JLFIN-l)  and 
TL( JLFIN) ,  if  all  the  excess  pore  pressures  have  absolute 
values  smaller  than  0.1  lb/ft2,  the  run  will  terminate. 

£.  KODPPR=C.  The  initial  excess  pore  pressures,  prior  to 
the  beginning  of  the  analysis,  are  constant;  i.e.,  inde¬ 
pendent  of  depth. 

£.  QTDATA=D.  If  the  run  is  terminated  because  all  pore  pres¬ 
sures  have  absolute  values  smaller  than  PPLIM,  then  de¬ 
tailed  data  will  be  output  at  the  final  time  of  analysis. 

_t.  K0D0RD=V.  The  consolidation  curve  is  defined  using  void 
ratios. 

£.  K0DINT=N.  The  linear  e-o  curve  is  assumed  between 
specific  points  that  define  this  curve. 

v.  K0PSP=C.  The  coefficient  of  consolidation  is  constant. 

Output 

41.  The  complete  output  file  is  shown  in  Table  21 .  Note  that 
values  less  than  zero  have  been  printed  in  the  file  of  echo  prints  where 
no  value  was  defined  originally;  e.g.,  when  input  of  a  set  of  data  is 
terminated  by  use  of  a  value  -1.0  for  tiie  appropriate  variable. 


Comparison  with  Hand  Solutions 


42.  With  this  simple  example,  a  hand  comparison  can  easily  be 
made  (Table  22).  Olson  and  Ladd  explore  the  comparison  of  classical  and 
the  finite  difference  analysis  in  more  detail. 


r 


Table  21 

Output  Data  File  for  Program  FD31  Example  Problem 


CVP.PKP 
*.05  .05 

EXAMPLE  PROBLEM  FOR  SETTLEMENT  PACKAGE 


TABLE  1 


ORIGIN 

A  L  LA 

V  E  R  CO 

N  D  I  T  I  0  N 

ORIGINAL 

NUMBER  OF 

SUBMERGED 

LAYER 

THICKNESS  NODES 

UNIT  WEIGHT 

NO. 

FEET 

(PSF  ) 

1 

1.00 

Z 

50.0 

2 

10.00 

11 

50.0 

TABLE  8 

DATA  ON  U 

A  T  E  R 

TABLE 

AND  D  R  A 

THE  UATER  TABLE  UILL  REMAIN  AT  A  CONSTANT  ELEVATION  WITH 
ZUTO  -  -30.0  FEET. 


BOTH  HORIZONTAL  BOUNDARIES  ARE  FREELV  DRAINING. 


THE  PROGRAM  HILL  TREAT  ANY  LAYER  AS  FREELY  DRAINING  IF 
CV.CE.  O.10E  03  SO. FT. PER  DAY. 


STATIC  PORE  UATER  PRESSURES  ABOVE  THE  UATER  TABLE  ARE  ASSUMED 
TO  BE  ZERO  AT  ALL  TIMES. 


ALL  SAND  LAYERS  ARE  ASSUMED  TO  DRAIN  FREELY  IN  THE  HORIZONTAL 
DIRECTION  AND  THUS  HAVE  ZERO  EXCESS  PORE  PRESSURES  AT  ALL  TIMES. 


TABLE  3 

INITIAL  EXCESS  PORE  PRESSURES 


LAYER 

NODE 

PORE  PRESSURE 

NO. 

NO. 

(PSF  ) 

1 

1 

0. 

i-a 

2 

0. 

2 

3 

0. 

2 

4 

0. 

2 

S 

0. 

2 

6 

0. 

2 

? 

0. 

2 

a 

0. 

2 

9 

0. 

2 

19 

0. 

2 

11 

0. 

2 

12 

0. 

(Continued) 
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Table  21  (Continued) 


TABLE 

4 

LOAD  HI 

S  T  0  R  V 

FILL 

TOTAL 

TINE  FILL 

TINE  FILL 

LAVER 

UNIT  WEIGHT 

CONSTRUCTION  BEGINS 

CONSTRUCTION  ENDS 

NUMBER 

(PCF ) 

(DAVS  ) 

( DAVS ) 

1 

112.4 

0. 

0. 

EMBANKMENT 

TINE 

HEIGHT 

(DAVS) 

(FT) 

0. 

0. 

0. 

20.00 

S000 . 

20.00 

TABLE  5 

TINES  F 

OR  OUTPUT 

0.1 

DAVS 

1.0 

DAVS 

4.0 

DAVS 

10.0 

DAVS 

40.0 

DAVS 

100.0 

DAVS 

400.0 

DAVS 

1000.0 

DAVS 

2000.0 

DAVS 

(Continued) 
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Table  21  (Continued) 


TABLE  6 


E  -  P  CURVES 
EFFECTIVE 


STRESS 

VOID 

LAVER 

(PSF) 

RATIO 

1 

0. 

1.0000 

10000.0 

1.0000 

2 

0. 

2. 0075 

2300.0 

1.9500 

SLOPE  OF 
REBOUND 
E  -  LOG  P 
LAVER  CURVE 

1  e.iE-ea 

2  0 . IE-02 


VALUES  OF  E  ARE  FOUND  BV  INTERPOLATING 
E  AND  P  NATURALLV. 


TABLE  7 


TABLE  OF  INPUT  VALUES  OF  CV  AND  PIC 


EFFECTIVE  COEFF.OF 
LAVER  STRESS  CONSOLIDATION 
NO.  PSF  SQ.FT/DAY 

1  NOT  INPUT  0.20E  03 

2  NOT  INPUT  0.S0E-01 


COEFF.OF 
PERMEABILITY 
FT/DAV 
0.20E  03 
0.50E-O1 


(Continued) 


Table  21  (Continued) 


TABLE  9 

TABLE  OF  CONTROL  PARAMETERS 


THE  MAXIMUM  UAIUES  OF  ALPHA  ARE  LIMITED  TO  O.S  DURING  LOADING  OR 
UNLOADING,  TO  B.5  FOR  ANV  NON-LOADING  PERIOD  EXCEPT  THE  LAST  ONE 
AND  TO  BETWEEN  0.S  AND  28.0  DURING  THE  LAST  LOADING  PERIOD. 

THE  GAUS-SEIDEL  ITERATIONS  (JILL  CONTINUE  UNTIL  NO  PORE  PRESSURE 
CHANGES  BV  MORE  THAN  0.10  PSF  FROM  ONE  ITERATION  TO  THE  NEXT  BUT 
THE  ANALYSIS  UILL  TERMINATE  IF  100  ITERATIONS  ARE  PERFORMED  FOR 
ANY  ONE  SET  OF  PORE  PRESSURE  CALCULATIONS. 

IF  THE  TM  TIME  USED  BETWEEN  ONE  OUTPUT  TIME  AND  THE  NEXT  EXCEEDS 
S.0  SECONDS,  THE  ANALYSIS  IS  TERMINATED  AFTER  OUTPUTTING  THE  DATA. 

IN  SUBROUTINE  CONSOL,  THE  PROGRAM  ITERATES  ON  THE  EFFECTIUE  STRESS 
EOUATION  UNTIL  NO  UALUE  OF  Pd)  CHANGES  BY  MORE  THAN  0.10  PSF. 

THE  NUMBER  OF  SUCH  ITERATIONS  IS  OUTPUT  AS  CHOP. 

IF  THE  MAXIMUM  FRACTIONAL  CHANGE  IN  ANY  UALUE  OF  CU  EXCEEDS  1.0 
(OR  IS  LESS  THAN  THE  RECIPROCOL  OF  THIS  NUMBER  FOR  DECREASING  VALUES 
OF  CU)  THEN  THE  PROGRAM  REDUCES  THE  TIME  STEP  TO 
DT*DT*9.9«CHCLIH/CHGMAX  AND  STARTS  ON  A  NEW  SET  OF  CALCULATIONS. 

FOR  THIS  PROBLEM  FCU-1.0  IF  THE  MAXIMUM  FRACTIONAL  CHANGE  IN  CU 
IS  LESS  THAN  0.1  THEN  THE  ANALYSIS  DOES  NOT  CYCLE  BACK.  IF  THE 
MAXIMUM  FRACTIONAL  CHANCE  IS  BETWEEN  0.1  AND  1.0  A  SET  OF  AVERAGE 
CU  VALUES  ARE  CALCULATED  AND  THE  CALCULATIONS  ARE  REPEATED. 

TABLE  10 


I  N  1  T 

I  A  L  C 

0  N  D  I  T 

IONS 

LAVER 

NODE 

2d) 

SIGNSOd) 

STATPPd  ) 

PREPPd  > 

Pd) 

S(I) 

NO. 

NO. 

PSF 

PSF 

PSF 

PSF 

PSF 

FT 

1 

1 

0. 

0. 

1872.00 

0. 

0. 

0. 

1-2 

2 

1.000 

112.40 

1934.40 

0. 

56.00 

0. 

2 

3 

2.000 

224.80 

1996.80 

0. 

100.00 

0. 

2 

4 

3.000 

337.20 

2059.20 

0. 

150.00 

0. 

2 

S 

4.000 

449.60 

2121.60 

0. 

200.00 

0. 

2 

e 

5.000 

S62.00 

2184.00 

0. 

250.00 

0. 

2 

7 

6.000 

674.40 

2246.40 

0. 

300.00 

0. 

2 

8 

7.000 

786.80 

2308.80 

0. 

350.00 

0. 

2 

9 

8.000 

899.20 

2371.20 

0. 

400 . 00 

0. 

2 

10 

9.000 

1011.60 

2433.60 

0. 

450.00 

0. 

2 

11 

10.000 

1124.00 

2496.00 

0. 

500.00 

0. 

2-3 

12 

11.000 

1236.40 

2558.40 

0. 

550.00 

0. 

LAVER 

> 

IN 

ZIN(IH) 

PIN< IN ) 

E0( IN ) 

CU( IN) 

PK< IN) 

DZ( IN  > 

DZO( IN 

NO. 

NO. 

FT. 

PSF 

SQ.FT/DAV 

FT'DAV 

FT. 

FT. 

1 

1 

0.500 

25.00 

1.0000 

8.260E  03 

0.206E  03 

1.690 

0.5000 

2 

2 

1.500 

75.00 

2.0056 

0.500E-01 

9 .5606-0 1 

1.000 

0.3327 

2 

3 

2.50O 

125.00 

2.0044 

0.500E-01 

0.S00E-01 

1.000 

0.3328 

2 

4 

3.500 

175.00 

2.0031 

0.S00E-91 

0. 5006-01 

1.000 

0.3330 

2 

5 

4.500 

225.00 

2.0019 

0.59OE-01 

0 . 500E-01 

1.000 

0.3331 

2 

6 

5.500 

275.00 

2.0006 

0.569E-01 

0.500E-01 

1.000 

0.3333 

2 

7 

6.500 

325.00 

1.9994 

0.506E-01 

0 . 50OE-01 

1.000 

0.3334 

2 

8 

7.500 

375.00 

1.9981 

0.500E-01 

0.509E-01 

1.000 

0.3335 

2 

9 

8.500 

425.00 

1.9969 

0.566E-01 

0.500E-01 

1.000 

0.3337 

2 

10 

9.500 

475.00 

1 . 9956 

0.506E-61 

0.500E-01 

1.000 

0.3338 

2 

11 

10.590 

525.00 

1.9944 

0.S00E-61 

9.500E-91 

1.000 

0.3340 

(Continued) 
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Table  21  (Continued) 


LAVER 

NO. 

1 

i-a 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


L.  IN 

1  1 
2  3 
2  3 
2  4 
2  S 
2  6 
2  7 
2  8 
2  9 
2  19 
2  11 


TAKE  11 


s 

OLIM 

ION  I  N  F  0  R  1 

RATIO 

N 

TINE  (DAYS)  * 

0. 

1 

TINE  STEP (DAYS)  • 

0.1 

LOAD  0  (PSF )  • 

2248.00 

ELEVATION 

OF  UATER  TABLE  (FT. 

>  •  30. 

0 

ELEVATION 

OF  TOP  OF  FILL  (FEET)  •  19.99 

NUMBER  OF 

ITERATIONS  • 

2 

UPPER  LIMIT  ON  ALPHA  • 

0.580 

MAXIMUM  DEVELOPED  ALPHA  • 

0.805 

SICFU  (PSF)  • 

624. 

5 

NO 

.  OF  CYCLES  THRU  COEFF  > 

1 

TN 

TIME  (SECONDS)  - 

0.014 

CHOP  (PSF) 

■ 

0.000 

MP 

• 

1 

TOTAL  TOTAL  EXCESS 

EFF. 

SETTLE 

NODE 

DEPTH 

STRESS  PP 

PP 

STRESS 

MENT 

NO. 

FT. 

PSF  PSF 

PSF 

PSF 

FT 

1 

0.088 

2872.5  1872.5 

0. 

1000.0 

0.00B 

2 

1.808 

2984.9  1934.9 

0. 

1050.0 

0.008 

3 

2.004 

3097.1  2992.1 

995.0 

105.0 

0.004 

■  4 

3.804 

3289.5  3859.5 

1000.0 

150.0 

0.004 

S 

4.004 

3321.9  3121.9 

1800.0 

200.8 

0.004 

6 

5.804 

3434.3  3184.3 

1000.0 

250.0 

0.084 

7 

6.004 

3546.7  3246.7 

1000.0 

300.0 

0.004 

8 

7.004 

3659.1  3389.1 

1000.0 

350.0 

0.084 

9 

8.004 

3771.5  3371.5 

1000.0 

400.0 

0.004 

10 

9.004 

3883.9  3433.9 

1000.0 

450.0 

0.004 

11 

10.004 

3996.3  3491.3 

995.0 

505.0 

0.004 

12 

11.000 

4108.4  2558.4 

0. 

1550.0 

0. 

EFF. 

CV  CVNEU 

DEPTH 

STRESS 

VOID  SQ.FT/  SQ.FT/ 

K 

PM 

FT. 

PSF 

RATIO  DAY  DAY 

FT/DAY 

ALPHA 

PSF 

8. SI 

1025.8 

i .oeatststt******** 

0.20E  03 

20.000 

1025.0 

1.51 

577.5 

1.993  8.0508  0.0500 

0.50E-01 

0.005 

577.5 

2.S8 

127. 5 

2.884  8.8508  8.8508 

0.50E-01 

8.005 

127.5 

3.50 

175.0 

2.083  0.0500  0.0500 

0.50E-01 

0.005 

175.0 

4.50 

225.0 

2.882  8.0508  0.0508 

0.50E-0! 

0.005 

225.0 

5.50 

275.0 

2.001  0.0500  0.0500 

0.50E-01 

0.005 

275.0 

6.58 

325.8 

1.999  0.0588  8.0500 

0.50E-01 

8.005 

325.0 

7.50 

375.0 

1.998  0.0500  0.0500 

0 . 50E-01 

0 . 805 

375.0 

8. SO 

425.8 

1.997  0.0500  0.0500 

0.50E-01 

0.085 

425.0 

9.58 

477.5 

1.996  8.0500  0.0500 

0.S0E-01 

8.885 

477.5 

18.50 

1027.5 

1.982  8.8500  0.0500 

0  50E-81 

8.005 

1027.5 

LAYER  COMPRESSION 
NO.  (INCHES) 

1  -0.0000 

2  8.1910 


( Continued) 


DZ 

FT. 

1.880 

8.9W 

1.888 

1.800 

1.880 

1.888 

1.088 

1.888 

1.088 

1.880 

8.996 
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Table  21  (Continued) 


TABLE  11 


SOLUTION  INFORMATION 


LAVER 

NO. 

1 

1-2 


TIME  (DA VS)  ■  1.® 

TINE  STEP l DAYS )  •  0.9 

LOAD  0  (PS F)  •  2248.00 

ELEVATION  OF  UATER  TABLE  (FT.)  •  30.0 

ELEVATION  OF  TOP  OF  FILL  (FEET)  •  19.99 
NUMBER  OF  ITERATIONS  •  3 


UPPER  LIMIT  ON  ALPHA  • 
MAXIMUM  DEVELOPED  ALPHA  * 
SIGFU  (PSF )  • 

NO.  OF  CYCLES  THRU  COEFF  • 
TH  TIME  (SECONDS)  • 

CHOP  (PSF)  ■ 

MP  - 


TOTAL 

TOTAL 

NODE 

DEPTH 

STRESS 

PP 

NO. 

FT. 

PSF 

PSF 

1 

0.009 

2872.6 

1872.6 

2 

1.009 

2985.0 

1935.0 

3 

2.005 

3097.1 

2949.2 

4 

3.005 

3209.5 

3058.3 

S 

4.005 

3321.9 

3121.9 

6 

5.005 

3434.3 

3184.3 

7 

6.005 

3546.7 

3246 . 7 

8 

7.005 

36S9.1 

3309.1 

9 

8.005 

3771.5 

3371.5 

10 

9.005 

3883.9 

3432.7 

11 

10.004 

3996.3 

3448.4 

12 

11.000 

4108.4 

2558.4 

0.500 

0.045 

624.6 

1 

0.271 

-0.000 

1 

EXCESS  EFF. 

SETTLE 

PP 

STRESS 

MENT 

PSF 

PSF 

FT 

0. 

1000.0 

0.009 

0. 

1050.0 

0.009 

952.1 

147.9 

0.O05 

998.8 

151.2 

0.005 

1000.0 

200.0 

0.005 

1000. 0 

250.0 

0.005 

1000.0 

300.0 

0.005 

1000. 0 

350.0 

0.005 

1000.0 

400.0 

0.005 

998.8 

451.2 

0.005 

952.1 

547.9 

0.004 

0. 

1550.0 

0. 

EFF. 

CV 

CVNEU 

DEPTH 

STRESS 

VOID 

SQ.FT/ 

SQ.FT/ 

K 

PM 

L. 

IN 

FT. 

PSF 

RATIO  DAY 

DAY 

FT/DAY 

ALPHA 

PSF 

1 

1 

0.51 

1025.0 

1.000***t*SttttttS* 

0.20E  03180.000 

1025.0 

2 

2 

t.st 

598.9 

1.993 

0.0500 

0.0500 

0.50E-01 

0.045 

598.9 

2 

3 

2.50 

149.5 

2.004 

0.0500 

0.0500 

0.50E-O1 

0.045 

149.5 

2 

4 

3.50 

175.6 

2.003 

0.0500 

0.0500 

0.50E-01 

0.045 

175.6 

2 

6 

4.50 

225.0 

2.002 

0.0500 

0.0500 

0.50C-01 

0.045 

225.0 

2 

6 

5.50 

275.0 

2.001 

0.0500 

0.0500 

0.50E-O1 

0.045 

275.0 

2 

7 

6.50 

325.0 

1.999 

0.0500 

0.0500 

0.5OE-01 

0.045 

325.0 

2 

8 

7.50 

375.0 

1.998 

0.0500 

0.0500 

0.50E-01 

0.045 

375.0 

2 

9 

8.50 

425.6 

1.997 

0.0500 

0.0500 

0.50E-O1 

0.045 

425.6 

2 

10 

9.50 

499.5 

1.995 

e.esee 

0.0500 

O.5OE-01 

0.045 

499.5 

2 

11 

10.50 

1048.9 

1.981 

0.0500 

0.0500 

0.50E-01 

0.045 

1048.9 

LAYER 

COMPRESSION 

NO. 

(INCHES) 

1 

0.0000 

2 

0.1098 

DZ 

FT. 

l.eee 

0.996 

1.000 

1.800 

1.008 

1.000 

1.000 

1.000 

1.000 

1.000 

0.996 


(Continued) 


Table  21  (Continued) 


TABLE  11 

SOLUTION  INFORMATION 

TINE  (DAYS)  •  10.0 

TINE  STEP (DAYS )  •  6.0 

LOAD  Q  (PSF )  •  2848.00 

ELEVATION  OF  WATER  TABLE  (FT.)  •  30.0 

ELEVATION  OF  TOP  OF  FILL  (FEET)  •  19.98 

NUNBER  OF  ITERATIONS  •  5 

UPPER  LIMIT  ON  ALPHA  -  0.516 

NAXINUN  DEVELOPED  ALPHA  •  0.303 

SICFU  (PSF)  •  624.9 

NO.  OF  CYCLES  THRU  COEFF  -  1 

TP1  TINE  (SECONDS)  •  0.013 

CHOP  (PSF)  •  0.000 

HP  -  1 


TOTAL  TOTAL  1 

EXCESS 

EFF. 

SETTLE- 

LAYER 

NODE 

DEPTH 

STRESS  PP 

PP 

STRESS 

NENT 

NO. 

NO. 

FT. 

PSF  PSF 

PSF 

PSF 

FT 

1 

1 

0.015 

2872.9  1872.9 

0. 

1000.0 

0.015 

i-e 

2 

1.015 

2985.3  1935.3 

0. 

1050.0 

0.015 

2 

3 

2.010 

3097.4  2665.7 

668.3 

431.7 

0.010 

2 

4 

3.008 

3209.7  2992.7 

933.0 

217.0 

0.008 

2 

S 

4.008 

3322.1  3111.9 

989.8 

210.2 

0.008 

2 

6 

5.008 

3434.5  3183.1 

998.6 

251.4 

0.008 

2 

7 

6.008 

3546.9  3246.5 

999.7 

300.3 

0.008 

2 

8 

7.00B 

3659.3  3307.9 

998.6 

351.4 

0.008 

2 

9 

8.008 

3771,7  3361.5 

989.8 

410.2 

0.008 

a 

10 

9.007 

3884.1  3367.0 

933.0 

517.0 

0.007 

2 

11 

10.006 

3996.3  3164.7 

668.3 

831.7 

0.006 

2 

12 

11.000 

4108.4  2558.4 

0. 

1550.0 

0. 

EFF. 

CU  CVNEU 

DEPTH 

STRESS 

VOID  SO. FT/  SQ.FT/ 

K 

pn 

DZ 

L 

.  IN 

FT. 

PSF 

RATIO  DAY  DAY 

FT/DAV 

ALPHA 

PSF 

FT. 

1 

1 

0.S2 

1025.0 

l.oietmmimm 

0.20E  03******I 

1025.0 

1.000 

2 

2 

1.51 

740.8 

1.989  0.0500  0.0500 

0.S0E-01 

0.303 

740.8 

0.994 

2 

3 

2.51 

324.3 

1.999  0.0500  0.0500 

0.50E-01 

0.300 

324.3 

0.998 

2 

4 

3.51 

213.6 

2.002  0.0500  0.0500 

0.S0E-01 

0.300 

213.6 

1.000 

2 

5 

4.51 

230.8 

2.002  0.0500  0.0500 

0.S0E-01 

0.300 

230.8 

1.000 

2 

6 

5.51 

275.9 

2.001  0.0500  0.0500 

0.50E-01 

0.300 

275.9 

1.000 

2 

7 

6.51 

325.9 

1.999  0.0500  0.0500 

0.50E-01 

0.300 

325.9 

1.000 

2 

8 

7. SI 

380.8 

1.998  0.0500  0.0500 

0.S0E-01 

0.300 

380.8 

1.000 

a 

9 

8.51 

463.6 

1.996  0.0500  0.0500 

0.50E-01 

0.300 

463.6 

1.000 

2 

10 

9.51 

674.3 

1.991  0.0500  0.0500 

0.50E-01 

0.300 

674.3 

0.998 

2 

11 

10.50 

1190.8 

1.978  0.0500  0.0500 

0.50E-01 

0.303 

1190.8 

0.994 

LAYER  COMPRESSION 
NO.  (INCHES) 


0.0000 

0.1821 


Table  21  (Continued) 


TABLE  11 


SOLUTION  INFORMATION 


TIHE  (DAVS)  • 

100.0 

TIHE  STEP (DAVS)  - 

a.o 

LOAD  0  (PSF )  • 

8248.00 

ELEVATION  OF  WATER  TABLE 

(FT.  ) 

•  30.0 

ELEVATION  OF  TOP  OF  FILL 

(FEET) 

•  19.96 

NUMBER  OF  ITERATIONS  • 

3 

UPPER  LIMIT  ON  ALPHA  - 

9.882 

MAXIMUM  DEVELOPED  ALPHA  • 

0.106 

SICFU  (PSF)  • 

626.7 

NO.  OF  CVCLES  THRU  COEFF 

• 

S 

TM  TIME  (SECONDS)  - 

0.06S 

CHOP  (PSF)  • 

0.906 

HP  •  1 


TOTAL 

TOTAL 

EXCESS 

EFF. 

SETTLE¬ 

LAVER 

NODE 

DEPTH 

STRESS 

PP 

PP 

STRESS 

MENT 

NO. 

NO. 

FT. 

PSF 

PSF 

PSF 

PSF 

FT 

1 

1 

0.043 

2874 . 7 

1874.7 

0. 

1000.0 

0.043 

1-2 

2 

1.043 

2987.1 

1937.1 

0. 

1050.0 

0.043 

2 

3 

2.03S 

3999.0 

2241.3 

242.3 

857.7 

0.035 

2 

4 

3.030 

3211.1 

2519.3 

458.2 

691.8 

0.030 

2 

5 

4.026 

3323.2 

2749.0 

625.7 

574.3 

0.026 

2 

6 

S.023 

3435.  S 

2916.2 

730.7 

519.3 

0.023 

2 

7 

6.021 

3547.7 

3014.1 

766.4 

533.6 

0.021 

2 

8 

7.019 

3660.0 

3040.7 

730.7 

619.3 

0.019 

2 

9 

8.017 

3772.2 

2997.9 

625.7 

774.3 

0.017 

2 

16 

9.913 

3884.4 

2892.5 

458.1 

991.9 

0.013 

2 

11 

10.607 

3996.5 

2738.7 

242.3 

1257.7 

0.007 

2 

12 

11.000 

4108.4 

2558.4 

0. 

1559.0 

0. 

L. 

IN 

depth 

FT. 

EFF. 

STRESS 

PSF 

CU 

UOID  SO. FT/ 
RATIO  DAV 

CVNEU 
SO. FT/ 
DAV 

K 

FT/DAV 

ALPHA 

PM 

PSF 

DZ 

FT. 

1 

1 

9.64 

1025.6 

1 .000** ************ 

0.20E  03394.248 

1025.0 

1.000 

2 

2 

1.54 

953.8 

1.984 

9.0500 

0.0506 

0.59E-01 

0.100 

953.8 

0.993 

a 

3 

2. 53 

774.7 

1.988 

0.9500 

0.0500 

0.S0E-01 

0.100 

774.7 

0.995 

a 

4 

3.53 

633.0 

1.992 

0.0500 

9.0500 

0.50E-01 

0.099 

633.0 

0.996 

a 

5 

4.52 

546.8 

1.994 

0.0500 

0.0500 

0.50E-01 

0.099 

546.8 

0.997 

a 

6 

5.52 

526.5 

1.994 

0.0500 

0.0500 

0.50E-01 

0.099 

526.5 

9.998 

a 

7 

6.52 

576.5 

1.993 

0.0500 

0.0500 

0.50E-01 

0.099 

576.5 

0.998 

a 

8 

7.S2 

696.8 

1.990 

0.0500 

0.0500 

0.50E-01 

0.099 

696.8 

0.997 

a 

9 

8. St 

883.1 

1.985 

9.0569 

9.0500 

0.S0E-01 

0.099 

883.1 

0.996 

a 

19 

9. SI 

1124.8 

1.979 

0.0500 

9.0509 

0.50E-01 

0.100 

1124.8 

0.995 

a 

11 

19.59 

1493.9 

1.972 

0.0500 

0.0S00 

0.50E-01 

0.100 

1403.9 

0.993 

LAVER  COMPRESSION 
NO.  (INCHES) 

1  9.0006 

2  9. 5180 


(Continued) 


Table  21  (Continued) 


TABLE  11 


SOLUTION  INFORMATION 


TIME  (DAVS)  • 

1000.0 

TIME  STEP (DAVS)  • 

22 

.0 

LOAD  0  (PSF )  • 

2248 

.00 

ELEVATION 

OF  UATER 

TABLE  (FT. 

)  •  30.0 

ELEVATION  OF  TOP  OF  FILL  (FEET)  *  19 

.92 

NUMBER  OF 

ITERATIONS  - 

3 

UPPER  LIMIT  ON  ALPHA  • 

4. 

314 

MAXIMUM  DEVELOPED 

ALPHA  • 

1. 

120 

SIQFU  (PSF)  • 

629.2 

NO.  OF  CVCLES  THRU 

COEFF  • 

11 

TM  TIME  (SECONDS) 

■ 

0. 

154 

CHOP  (PSF) 

• 

-0. 

000 

HP  • 

1 

TOTAL 

TOTAL 

EXCESS 

EFF. 

SETTLE¬ 

LAVER 

NODE  DEPTH 

STRESS 

PP 

PP 

STRESS 

MENT 

NO. 

NO.  FT. 

PSF 

PSF 

PSF 

PSF 

FT 

1 

1  B.B83 

2877.2 

1877.2 

0. 

1000.0 

0.083 

1-2 

2  1.083 

2989.6 

1939.6 

0. 

1050.0 

0.083 

2 

3  2.075 

3101.5 

2004.2 

2.8 

1097.2 

0.075 

2 

4  3. OSS 

3213.3 

2068.6 

5.3 

1144.7 

0.066 

a 

5  4.058 

3325.2 

2132.4 

7.2 

1192.8 

0.058 

2 

S  5.050 

3437.1 

2195.6 

8.4 

1241.6 

0.050 

2 

7  6.041 

3549.0 

2257.8 

8.9 

1291.1 

0.041 

2 

8  7.033 

3660.9 

2319.3 

8.4 

1341.6 

0.033 

2 

9  8.025 

3772.8 

2379.9 

7.1 

1392.9 

0.025 

2 

10  9.017 

3884.6 

2439.8 

5.2 

1444.8 

0.017 

2 

11  10.008 

3996.5 

2499 . 2 

2.7 

1497.3 

0.008 

2 

12  11.000 

4108.4 

2558.4 

0. 

1550.0 

0. 

L. 

IN 

DEPTH 

FT. 

EFF. 

STRESS 

PSF 

CV 

VOID  SQ.FT/ 
RATIO  DAV 

CVNEU 

SQ.FT/ 

DAV 

K 

FT/DAV 

ALPHA 

PM 

PSF 

DZ 

FT. 

1 

1 

0.58 

1025.0 

i.eetimmmim 

0.20E  03******* 

1025.0 

1.000 

2 

2 

1.58 

1073.6 

1.981 

0.0500 

0.0500 

0.50E-01 

1.120 

1073.6 

0.992 

2 

3 

2.57 

1121.0 

1.979 

0.0500 

0.0500 

o.soE-ei 

1.120 

1121.0 

0.992 

2 

4 

3.56 

1168.8 

1.978 

0.0500 

0.0500 

e.soE-oi 

1.120 

1168.8 

0.992 

2 

5 

4.55 

1217.2 

1.977 

0.0500 

0.0500 

0.50E-01 

1.120 

1217.2 

0.992 

2 

6 

5.55 

1266.3 

1.976 

o.osee 

0 . 0500 

o.soE-ei 

1.120 

1266.3 

0.992 

2 

7 

6.54 

1316.4 

1.975 

0.0500 

0.0500 

0.50E-01 

1.120 

1316.4 

0.992 

2 

8 

7.53 

1367.2 

1.973 

0.0500 

0.0500 

0.50E-01 

1.120 

1367.2 

0.992 

2 

9 

8.52 

1418.9 

1.972 

0.0500 

0.0500 

o.soE-ei 

1.120 

1418.9 

0.992 

2 

10 

9.51 

1471.1 

1.971 

0.0500 

0.0500 

0.50E-01 

1.120 

1471.1 

0.992 

2 

11 

10.50 

1523.6 

1.969 

0.0500 

0.0500 

0.50E-01 

1.120 

1523.6 

0.992 

LAVER  COMPRESSION 
NO.  (INCHES) 

1  0.0600 

2  0.9944 
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Table  21  (Continued) 


TABLE  11 


SOLUT 

I  0  N 

INFORM 

A  T  I 

0  N 

TINE  (DAVS)  ■ 

2000.0 

TINE  STEP (DAVS)  • 

90 

.4 

LOAD  0  (PSF )  • 

2248 

.00 

ELEVATION 

OF  UATER 

TABLE  (FT.) 

•  30.0 

ELEVATION 

OF  TOP  OF  FILL  (FEET) 

•  19 

.92 

NUMBER  OF 

ITERATIONS  ■ 

1 

UPPER  LIMIT  ON  ALPHA  ■ 

7. 

947 

MAXIMUM  DEVELOPED 

ALPHA  • 

4. 

598 

SIGFU  (PSF)  ■ 

629.2 

NO.  OF  CVCLES  THRU 

COEFF  • 

9 

TM  TIME  (SECONDS) 

• 

6. 

373 

CHOP  (PSF) 

■ 

-0. 

000 

HP  • 

1 

TOTAL 

TOTAL  EXCESS 

EFF. 

SETTLE¬ 

LAVER 

NODE  DEPTH 

STRESS 

PP 

pp 

STRESS 

MENT 

NO. 

NO.  FT. 

PSF 

PSF 

PSF 

PSF 

FT 

1 

1  •.•83 

2877.2 

1877.2 

0. 

1000. 0 

•  .083 

1-2 

2  1.083 

2989.6 

1939.6 

0. 

1050.0 

0.083 

2 

3  2. *75 

3101.5 

2001 . 6 

0.1 

1099.9 

0.075 

2 

4  3.067 

3213.4 

2063.6 

0.3 

1149.7 

•  .•67 

2 

5  4.058 

3325.2 

2125.6 

0.4 

1199.6 

•  .058 

2 

6  5. OS* 

3437.1 

2187.5 

•  .4 

1249.6 

•  -•50 

2 

7  6.042 

3S49.0 

2249.4 

•  .4 

1299.6 

•  .•42 

2 

8  7.033 

3660.9 

2311.3 

0.4 

1349.6 

•  .033 

2 

9  8.025 

3772.8 

2373.1 

0.3 

1399.7 

•  .•25 

2 

10  9.017 

3884.6 

2434.9 

0.2 

1449.8 

•  .•17 

2 

11  10.008 

3996.5 

2496.6 

0.1 

1499.9 

0.008 

2 

12  U.MO 

4198.4 

2558.4 

0. 

1559.0 

». 

EFF. 

CV 

CVNEU 

DEPTH 

STRESS 

VOID  SQ.FT/ 

SQ.FT/ 

K 

PN 

IN  FT. 

PSF 

RATIO  DAV 

DAV 

FT/D AV  ALPHA 

PSI 

1 

1 

0.58 

1025.0 

i.0eo»*t*s*t*t**s* 

0.20E  03******* 

1025.0 

2 

2 

1.58 

1*74.9 

1.981 

0.0500 

0.0500 

0.50E-01 

4.597 

1074.9 

2 

3 

2.57 

1124.8 

1.979 

0.0500 

0.0500 

0.50E-01 

4.597 

1124.8 

2 

4 

3.56 

1174.7 

1.978 

0.0500 

0.0500 

0.50E-01 

4.597 

1174.7 

2 

5 

4.55 

1224.6 

1.977 

0.0500 

e.esee 

0.50E-01 

4.597 

1224.6 

2 

6 

5.55 

1274.6 

1.976 

0.0500 

0.0500 

e.soE-ei 

4.597 

1274.6 

2 

7 

6.54 

1324.6 

1.974 

0.0500 

0.0500 

O.5OE-01 

4.598 

1324.6 

2 

8 

7.53 

1374.7 

1.973 

0.0500 

0.0500 

0.50E-01 

4.598 

1374.7 

2 

9 

8.52 

1424.7 

1.972 

0.0500 

0.0500 

0.50E-01 

4.598 

1424.7 

2 

10 

9.51 

1474.8 

1.971 

0.0500 

0.0500 

0.50E-01 

4.598 

1474.8 

2 

11 

10.5* 

1524.9 

1.969 

0.0500 

0.0500 

e.soE-ei 

4.598 

1524.9 

LAVER  COMPRESSION 
NO.  (INCHES) 

1  •.••09 

2  0.9997 


DZ 

FT. 

1.000 

•  .992 

•  .992 
0.992 
0.992 

•  .992 

•  ■992 
0.992 

•  .992 
0.992 

•  .992 


(Cont inued) 


Table  21  (Concluded) 


r 


TABLE  12 


R  V  0  F 

T  I  A  E  S  E 

TTIEBEN 

DECREE  OF 

TIRE 

SETTLEMENT 

CONSOLIDATION 

(DATS) 

(FEET) 

(PERCENT) 

B.  1 

0.00842 

10.1 

l.B 

0.00915 

11.0 

4.0 

0.01140 

13.7 

10.0 

0.01517 

18.2 

40.0 

0.02757 

33.1 

100.0 

0.04267 

51.2 

400.0 

0.07417 

89.0 

1000.0 

0.08287 

99.5 

2000.0 

0.08331 

100.0 

I 


85 


Table  22 

Comparison  of  FD31  and  Hand  Solutions  for 
Settlement  and  Degree  of  Consolidation 


Settlement, 

ft 

FD31  Solution 

0.08331 

Hand  Solution 

0.083 

Time 

Degree 

of  Consolidation 
percent 

days 

FD31  Solution 

Hand  Solution 

0.1 

10.1 

3 

1.0 

11.0 

5 

4.0 

13.7 

10 

10.0 

18.2 

16 

40.0 

33.1 

32 

100.0 

51.2 

51 

400.0 

84.0 

88 

1000.0 

99.5 

99 

2000.0 

100.0 

100 

86 
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APPENDIX  A:  INPUT  FOR  EMBANKMENT  LOADS— PROGRAM  10016 


1.  This  appendix  describes  the  additional  input  data  for  program 
10016  that  are  required  to  analyze  the  vertical  stress  in  the  foundation 
beneath  an  embankment.  Following  the  lines  describing  the  loaded  rec¬ 
tangular  ares  (if  any)  is  a  set  of  lines  describing  the  shape  and  weight 
of  the  embankment  loading  (see  Figures  Al  and  A2) .  These  lines  will  be 
necessary  only  under  the  option  where  KODE  is  input  as  2  or  3.  The  first 
line  required  to  describe  the  embankment  consists  of  the  following  input 
variables  (see  line  430  in  Table  Al) : 

a.  NCOR.  NCOR  is  the  number  of  pairs  of  X  and  Y  coordinates 
used  to  describe  the  cross  section  of  the  embankment  load¬ 
ing.  The  maximum  allowable  value  of  NCOR  is  25. 

_b.  GAMMA.  The  variable  GAMMA  is  the  unit  weight  of  the  embank¬ 
ment  fill  in  units  of  weight  and  length  compatible  with 
the  other  input  data. 

£.  THICK.  THICK  is  the  input  variable  which  determines  the 
number  of  layers  used  to  approximate  the  embankment  load¬ 
ings.  THICK  represents  the  maximum  allowable  thickness  of 
any  layer  used  in  the  approximation  of  the  embankment 
loading. 

d.  YMAX.  The  value  of  YMAX  is  the  longitudinal  distance  from 
the  cross  section  to  the  end  of  the  embankment  in  the 
positive  Y  direction.  YMAX  should  in  all  cases  be  greater 
than  or  equal  to  zero,  since  the  cross  section  (X-Z  plane) 
defining  the  embankment  loading  is  assumed  to  be  at  Y  =  0. 

e.  YMIN.  The  value  of  YMIN  is  the  longitudinal  distance  from 
the  cross  section  to  the  end  of  the  embankment  in  the  nega¬ 
tive  Y  direction.  YMIN  should  in  all  cases  be  less  than  or 
equal  to  zero  since  the  cross  section  (X-Z  plane)  defining 
the  embankment  loading  is  assumed  to  be  at  Y  =  0. 

2.  The  remaining  cards  required  to  describe  the  embankment  loading 
consist  of  a  series  of  lines  each  defining  a  pair  of  corner  points  (X,Z) 
of  the  embankment  cross  section.  The  number  of  these  lines  will  corre¬ 
spond  to  the  value  of  NCOR  described  above.  These  corner  point  lines 
should  be  input  in  the  same  sequence  as  they  appear  in  the  embankment 
cross  section  (see  Figure  A2) .  In  other  words,  the  input  sequence  should 
be  the  same  as  would  be  found  by  proceeding  around  the  perimeter  of  the 
embankment  cross  section  in  either  a  clockwise  or  a  counterclockwise 


Figure  A2.  Cross  section  and  plan  view  of  embankment  for  input  description 
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Table  A1  (Concluded) 


DO  YOU  WISH  TO  RUN  PROGRAM  FROM  EXISTING  DATA  FILE? 

=ye: 
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DO  YOU  WANT  OUTPUT  WRITTEN  TO  AN  OUTPUT  FILE? 
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CONFIGURATION. 
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direction.  The  lines  required  for  each  corner  point  will  have  a  format 
described  as  follows  (see  lines  440-470  in  Table  Al) : 

a.  X(I) .  X(I)  is  the  X  coordinate  of  a  corner  (break  point) 
in  the  shape  of  the  embankment  cross  section. 

b.  Z(l) ■  Z(I)  is  the  Z  coordinate  of  a  corner  (break  point) 
in  the  shape  of  the  embankment  cross  section.  The  Z  coor¬ 
dinates  are  referenced  to  the  lowest  Z  coordinate  which 
has  a  value  of  zero. 
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